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sense that whenever such scattering events occur, nonperturbative quantum few-fermion scattering 
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properties of the wave function are studied, and its parameters in the Bose-Einstein condensed limit 
are computed beyond mean-field. These results enable us to predict many observable properties of 
this Fermi gas with well-controlled accuracy, such as the superfiuid pairing function, the four-fermion 
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I. INTRODUCTION 

In this paper we consider a superfluid Fermi gas with 
strong enough attractive interactions such that molecules 
of pairs of fermions are formed and Bose-Einstein con- 
densed. A physical example of the system we have in 
mind is a two-component Fermi gas near a Feshbach res- 
onance, with positive inter-fermionic scattering length. 
Our goal is to develop a new form of perturbation theory 
which will enable us to address fundamental and, now, 
experimentally accessible properties of this gas. One can 
not determine many observables of this system by sim- 
ply resorting to the common wisdom of conventional Bose 
gases. In many instances, the underlying fermionic char- 
acter must play a profound role. The starting point for 
the present paper is a wavefunction ansatz introduced 
earlier [3| for the ground state of a superfluid Fermi gas. 
In this earlier work we explored the equation of state 
and its relation to breathing mode experiments. Here 
we present in considerable detail our theoretical method- 
olgy. To demonstrate the predictive power of our the- 
ory, we also compute the superfluid pairing function and 
the four-fermion correlation function. In future work we 
will address the fermionic momentum distribution, the 
two-body reduced density matrix, the higher correlation 
functions as well as other properties. 

Because the attractive interactions are strong enough 
to form bound states, whenever two or more molecules 
collide, the few-body physics is highly nonperturbative. 
In this paper we build on the observation that there is, 
nevertheless, an underlying perturbative component: the 
interactions are weak in the sense that the probability 
amplitudes (or component terms of the many-body wave 
function) associated with collisions are small, provided 
that the density of the system is relatively low. This 



is precisely the situation in which standard many-body 
perturbation techniques, based on weak inter-particle in- 
teractions are inappropriate. By contrast, here we de- 
velop a new kind of perturbation method, taking full ad- 
vantage of these small amplitudes. Whenever molecules 
approach each other, however, the behavior of the con- 
stituent particles is treated nonperturbatively. In this 
way a complete many-body theory fully compatible with 
the nonperturbative few-body physics will be presented. 

Closely related ideas have been used in the companion 
paper '5| to address the problem of a gas of point-like par- 
ticles. This is, in some sense, the simpler analog of the 
bound states problem. In this companion paper |2| the 
ground state of the dilute Bose gas of structureless parti- 
cles with quite arbitrary interactions is studied, starting 
from a general many-body wave function: 

|^boson)=exp(£l ^ <...k,4.---<)|0>, 

>=l^'kl...kp ^ 

where the &'s are boson annihilation operators, the k's 
are momenta, and |0) is the particle vacuum. To our sur- 
prise, this familiar system has some fundamental proper- 
ties unknown to people, because the prevailing theories to 
date (Bogoliubov theory, pseudopotential method, effec- 
tive field theory, etc) are all low-energy effective theories, 
capable primarily of describing the physics at the length 
scale of the superfluid healing length. In the companion 
paper such a limitation is removed, and the short- 
distance structure is clarified. As a byproduct, some fun- 
damental theorems derived initially for a Fermi system 
^ 01 are successfully extended to a dilute Bose system 

In this paper we apply the same approach to a Bose- 
Einstein condensate of molecules formed out of strongly 
attractive fermions. Each molecule is a bound state of 
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two fermions in two spin states (or any two internal 
states) t= +1 and J,= — 1. In Ref. we presented a 
generalized many-body wavefunction of the form 

IV'fermion) = SXp | ^ ^ O^^Ik^ • • 4i • ■ • 4^ |0> > 

_ (1) 

where p = 2, 4, 6, • • • must be even [3, the subscript K.^ 
is the shorthand for kitJi, and a labels the spin state. 
The coefficient Q;j^^...j^ has some elementary properties: 
1) it is antisymmetric under the interchange of any two 
subscripts, 2) it is zero whenever the sum of all its sub- 
scripts is nonzero, and 3) it is zero whenever the sum of 
a nontrivial subset f Appendix of its subscripts is zero 
(see Sec HTEl) . (Here the sum of subscripts is defined by 
summing the momenta and the associated cr's separately, 
and when both sums are zero, the sum of these subscripts 
is called zero.) 

The well-known BCS or Eagles-Leggett (EL) wave 
function of fermionic superfluids is the lowest order ap- 
proximation of Eq. . This wavefunction and its impli- 
cations for experiment have been reviewed rather exten- 
sively elsewhere @. Retaining only the p = 2 term, one 
can verify that 1-0) becomes equivalent to a more familiar 
form of the EL wave function: 

IV'El) ^("k + ^'k4Tciki)|0)■ 
k 

In Sees, ini mil and IIVI the ideas and methods in 
Ref. ^ are extended to bosonic molecules. Here we 
present details to support our earlier conclusions for 
the ground state energy density calculated to the 2.5-th 
order in the molecular density. As discussed in Ref. 0, 
the result indicates that the Lee- Yang formula V] is valid 
not only for structureless bosons, but also for composite 
ones. 

In Sec.|n]we define the system, discuss those properties 
of the many-body wave function that are valid for all 
densities, and introduce a new diagrammatic method. 

In Sec. mil we discuss the few-body functions and the 
properties that will be referred to in the subsequent 
many-body low-density expansions. 

In Sec. IIVI we introduce our low-density expansion 
method, derive the power-counting formula for such ex- 
pansions, and compute the parameters of the wave func- 
tion beyond mean-field. The equation of state is also de- 
termined beyond mean-field. Our method is variational, 
but unlike the traditional sense of this term, we identify 
all the significant contributions up to any given order in 
the density, ensuring that our results are asymptotically 
exacts order by order. 

In Sec. |Vlwe compute two other physical observables, 
using the established knowledge of the wave function, and 
discuss the motivation and the feasibility of computing 
some other fundamental observables. 

In Sec. IVII we outline the extension of the calculation 
to third order in the molecular density, and predict that 



the Wu term 8] must also be present, and contain the 
same scattering length as the Lee- Yang term 0. 

With each increasing order in the density expansion of 
the energy, the many-body wave function itself is deter- 
mined more and more accurately. In this way, we can 
determine all other physical observables with a high de- 
gree of accuracy by evaluating their expectation values 
within the many-body wave function. This approach 
also makes it possible to address quantitative properties 
slightly above zero temperature. 

Some trends in the system, upon increasing density, are 
clearly visible in the low-density expansion. This opens 
a new avenue towards the many-body physics in the uni- 
tary scattering regime. 



II. FERMIONIC SUPERFLUID: SOME 
GENERAL ANALYSES 

All results in this section are completely independent of 
the density, and equally valid in regimes such as the uni- 
tarity limit of a Fermi gas with large scattering length. 
They are however a prerequisite of the low-density ex- 
pansion in Sec. IIVI In light of the possibility of more 
concrete progress in the high density regime, it is de- 
sirable to separate the generally applicable results from 
those that are restricted to low density. 



A. Hamiltonian 

We study a system of fermions with identical mass m 
(let h = m = 1) in two internal states (hereafter simply 
called spin states). For simplicity we restrict ourselves 
to a system in which there is only interaction between 
fermions in the opposite spin states (t and |). The hamil- 
tonian H and the operator N for the total number of 
fermions satisfy 

H-fiN^Y. f dV4(r)(-VV2-M)c,(r) + 

y"d^ri • • •d^r4C/(rir2r3r4)cj(ri)cJ(r2)c|(r4)c|(r3), 

(2) 

where /i is the fermionic chemical potential, a =t, i, and 
Ccr(r) is the standard fermion annihilation operator in the 
coordinate representation, satisfying {co.(r), c^, (r')} = 

We restrict ourselves to an interaction U with the fol- 
lowing properties: 

• translational symmetry, U{ri + r, r2 + r, + r, r4 + 
r) = C/(rir2r3r4), for any r. 

• rotational symmetry, U{Tlri,TZr2,TZr3,TZr4) = 
C/(rir2r3r4), for any rotational transformation TZ. 
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Galilean symmetry, t/(rir2r3r4) 
r2 ^ ra + r4. 



whenever ri - 



• finite-rangedness: there is a finite length scale I 
such that C/(rir2r3r4) approaches zero sufficiently 
fast whenever any \ri — rj\ ^ I, where 1 < «, j < 4. 

• the interaction supports a nondegenerate s-wave 
bound state (the molecule) of two fermions, whose 
binding energy in the absence of any other parti- 
cles will be denoted by Em < (independent of the 
center-of-mass velocity of the molecule, due to the 
Galilean symmetry); the typical distance between 
the two fermions in the molecule will be denoted 
by Tm', the zero- velocity scattering length between 
two molecules, a„i, is positive and finite. 

• lower bound states do not exist, or they do exist 
but their formation rates are so low that we can 
omit them Q. 

Also, H ^ W ^ L/(rir2r3r4)* = [/(r3r4rir2). 

For later convenience, we use a symbol R = (ra), 
called "spin-spatial vector"; and similarly K = (ktr), 
called "spin-momentum " . We also define some opera- 
tions involving them. The "norm-square": = k^. The 
"addition" : Ki +K2 = (ki -|-k2, cti -1-0-2), where each t is 
counted as +1 and each J, is counted as —1; for any even 
number of spin-momenta Ki • • • Kp, if both ai + - ■ ■+(Tp = 
and ki -I- • • • -I- kp = 0, we say that Ki + • • • + Kp ~ 
0. The sum of a spin-momentum and a momentum: 
K -I- q = (k -|- q, (t). The "opposite spin-momentum": 
-K = (-k, -a). The "integration": / d^i? = Y.a ! d^^' 
and j6?K/{2T:f = / d^fc/(27r)^ and similarly a 
summation: = X^o-k ~ ^/ d^i4r/(27r)'^, where Vl is 

the system's volume. Finally, if both cti + • • • + cTp = 
and ki -)- • • • -1- kp ~ g, we say that Ki + • ■ ■ + Kp ^ q, 
where q is some small momentum scale. 

Now we define the function [/(R1R2R3R4), sat- 
isfying: J7(RiR2R3R4) = -[/(R2R1R3R4) = 
— C/(RiR2R4R3); it is zero whenever cti + (72 7^ or 
CT3 + cr4 0; U{vi t,r2 i,r3 t,r4 i) = U{viV2V:iYi). 
Equation Q is thus rewritten as 

H-^iN= j d3i?c^(R)(-VV2- At)c(R) 

+ i y A^Ri ■ ■ ■ d3i?4J/(RiR2R3R.4) 
X c^(Ri)c^(R2)c(R4)c(R3). 
Fourier transforms: 

Ck^ = n-^/"^ j dV c(rCT) exp(-ik • r). 



TJ — 
*-^kl(7lk2CT2k3(73k4CJ4 



d^ri • • • d^r4 



X [/(riCrir2Cr2r30-3l"4CT4) 

X exp(-iki • ri - ik2 • V2 + iks • + ik4 • r4). 



K 

+ iXl^KiKaKakiCKiCic^CK.CKg. (3) 
K's 

{ck'Ck^} = ^KK' and {cj^.c^,} = 0. 

C^KiK2K3K4 = unless Ki + K2 = K3 + K4. When 
this spin-momentum conservation condition is satisfied, 
[/ is a smooth, function of the remaining three indepen- 
dent momenta, as a consequence of the condition that 
the interaction is finite-ranged. Third, 



K1K2K3K. 



= -u, 



K2K1K3K4 



K1K2K4K3 



f^K3K4KiK2 ~ C^Ki+k,K2+k,K3+k,K4-l-k, 



where the last equality is due to the Galilean symmetry 
(k is an arbitrary momentum). Fourth, in the thermo- 
dynamic limit, any nonzero C/K1K2K3K4 scales like [l0| 



K1K2K3K4 



In this paper, an expression like x ~ f2'5 only refers to 
the behavior of xiorVL^Q!] since the often dimensionful 
coefficient (independent of SI) is suppressed, one should 
never use such expression to infer the dimension of x. 



B. Wave function 

The ground state wave function we have adopted is 

|V') = exp|^^^a(f)^|0), 
where p = 2, 4, 6, • • • must be even (5|, and 

^' Ki - Kp 

is called the (p/2)-th order generator; 



' Ki-Kp 

is called the (p/2)-th order degenerator. The prefactor 
1/pl is for the antisymmetry of the c^^'s (or c's). 

In a very similar manner as in j2| , we can establish the 
following properties of ct^\...j^^- 

Firstly, we assume that the ground state is translation- 
ally invariant, and that the superfluid pairing occurs in 

(v) 

the spin-singlet channel only. So Q^k ...k ~ ^ whenever 
Ki + ---+Kp^O. 

Secondly, in the thermodynamic limit, in which the 
density of fermions 

n = N/n 
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is kept constant, any nonzero coefficient cx^[...j^ scales 
with the system's volume as 

^l.-K.-f^^-"/^- (4) 

(p) 

Thirdly, c^k • k antisymmetric under the inter- 
change of Ki and K^ , where I < i < j < p. 

C. Diagrams 

This subsection is largely parallel to a similar subsec- 
tion in 0. The only essential additional feature is due 
to the antisymmetry of the fermionic wave function. 

The evaluation of any physical observables can be re- 
duced to the following general problem. 

Consider an arbitrary product O = 01O2 ■ ■ ■ ol, where 
each Oj (1 ^ i ^ L) is a basic fermionic operator (cre- 
ation or annihilation) in the spin-momentum space. The 
ground state expectation value 




is calculated using the Wick's theorem. We assign each 
fermion operator in the above expression a "time" : the 
creation operators in the generators are assigned an ear- 
liest time ig, the annihilation operators in the degen- 
erators a latest time tj_, and Oi through are as- 
signed some intermediate times, ti through t^, such that 
td > ti > ■ ■ ■ > II > tg. The actual magnitudes of these 
"times" do not matter; only the "causal order" matters. 
Equation (|SJ) is thus rewritten as 

{0\TO expfer=2 + E;?=2 "^^'') |0> 
(0) = ; (6) 

(o|rexp(Er=2"(^^^ + E;^=2«(^'))|o) 

where T is the "time" ordering operator. 

The basic formulas to be used in the evaluation of 
Eq. © are 

{0\Tc^tCK't'm = mclA't' |0> = 0, (7a) 



(0|TcKt4,,,|0) = -{0\Tcl,,,c^,\0) = SKK,0it - t'). 

(7b) 

Equation © is then expressed as an infinite sum of 
all the topologically distinct diagrams with the following 
rules. 

• Each Oi is represented a small solid circle if it is 
a fermion creation operator, or a small hollow one 
if it is an annihilation operator, linked by a single 
solid line, with the designated spin-momentum. 



• Each p/2-th order generator (degenerator) is rep- 
resented by a large solid (hollow) circle, linked by 
p solid lines; it contributes a factor (or its 
complex conjugate). 

• Each solid line is associated with a single spin- 
momentum, and it can only connect an earlier solid 
object with a later hollow one, in accordance with 
Eq. The diagram is thus bipartite. 

• Impose spin-momentum conservation on all the 
internal vertices, that is, the sum of all 
the spin-momenta leaving/entering any genera- 
tor/degenerator must equal zero. Sum over all the 
independent loop spin-momenta. 

• Each diagram has an overall symmetry factor 1/5, 
and S is the number of ways of interchanging com- 
ponents without changing the diagram. 

• Vacuum bubbles (disconnected parts which are not 
connected to any external vertex) are all canceled 
by the denominator in Eq. ©. 

The fermionic operators anticommute under the 
"time" -ordering operator; the subscripts of Q!k|...k 
anticommute. It is easy to show that there is a simple 
method to determine the sign of a diagram. 

Step 1. Divide all the solid lines of the diagram into 
a collection of continuous trajectories. Each solid line 
can only be contained by one trajectory, and it can ap- 
pear in this trajectory only once. Each trajectory is ei- 
ther an external trajectory, or an internal one. Each 
external trajectory starts from an external vertex o'*, 
ends at a different external vertex o'^, and is of the form 
i;oKiUiK2f;2 • • ■ Km^m, where M > 1, uo = o^, vm = o*", 
and the remaining tj's are internal vertices. K's are spin- 
momenta, that is, the solid lines. Each internal trajectory 
starts from an internal vertex vq , and must come back to 
this same vertex; it is of the form woKi'i;iK2W2 • • • ICmvm, 
where vm = vo, and M is necessarily a positive even 
number (recall that the diagram is bipartite); an internal 
trajectory can never pass an external vertex. The total 
number of external trajectories is necessarily L/2. The 
L external vertices is thus divided into L/2 pairs, and 
can be written as o\o\o'!2p\ ■ ■ ■ 0^/2*^2/2' suppose that it 
takes P interchanges of pairs of fermionic operators to 
transform the sequence oi • • • ol to this new one. Sup- 
pose that X of the L/2 external end points (o'^'s) are 
fermion annihilation operators, that is, hollow circles. 

Step 2. For every p-line internal vertex u, write a sepa- 
rate symbol a^^^ if it is solid, or a^^)* if it is hollow; leave 
p blank positions for its subscripts. 

Step 3. To each internal vertex v passed through by 
each trajectory, assign a pair of ordered subscripts KK', 
where the 3-tuple KwK' (exactly in this order) is a seg- 
ment of the trajectory. If v is the starting-ending point of 
an internal trajectory (containing M solid lines), assign 
to it a pair of ordered subscripts KiKm- Multiply the 
final expression by an overall factor (— 1)^"*"^. 
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FIG. 1: A diagram in the expansion of {cj^-^Cj^^Cj^^Cj^^) ■ 

There are different ways to divide the same diagram 
into trajectories, but the resuhs are equal. 

We iUustrate these rules with a term in the expan- 
sion of {cIq^cI^^cI^^Ciq^) , shown in Fig. ^ The diagram 
can be divided into two trajectories: 01K.1AK.5BK.4Oi 
and 02K2AK6BK7CK303, so (-1)^ = +1, X = 1, 
and (— — —1. The expression with the cor- 
rect sign is therefore -AKiKsKsKe^KsKiKsKTCKTKs = 
^KiKsKsKs^KiKyKsKeCKaKr- Incorporating the other 
rules, we get 

+ 2 X] ^KiK2K5K6/5K4,-K3,K5K6'^k3 7 

where we have given special names to the coefficients of 
the first three generators: 

"K = "k!-K' (8a) 

/^KiKsKaKi = Q^KlKaKaKi' (8^) 

— (6) /n \ 

7K1K2K3K4K5K6 = "KiKaKaKiKsKe,- l^^j 

There are quite a few other ways to divide 
Fig. n] into trajectories. For instance, 01K1AK2O2, 
o^K^CKjBKiOi, and AK^BK^A, where the last tra- 
jectory is an internal one; so (—1)^ = +1, X = 1, and 
(— 1)"^+^ = —1. The expression with the correct sign 
is therefore -AKiKjKsKgSKrKiKsKeC'KaKy, equivalent 
with the one we wrote above. Even if we reverse the di- 
rection of any trajectory, the final result is still the same. 

D. Diagrams for E — fiN 

Like in 0, we use a half-solid-half-hollow diamond 
to facilitate the diagrammatic calculation of i? — fxN — 
{H — /iiV). If each side is attached by a single solid line, 
the diamond contributes a factor — /i), where K 

is the spin- momentum flowing through it. If each side 
is attached by two solid lines, the diamond contributes 
a factor ?7kiK2K3K4, where K1K2 leave the diamond's 
solid side, and K3K4 enter its hollow side. 

The above rules are very similar to those in We 
now describe how to determine the sign of each diagram. 

The solid lines in such a diagram can still be divided 
into a collection of trajectories. 

For a trajectory vqKiVi ■ • ■ Kmvai (vq = vm) which 
does not involve the diamond, the rule is the same as in 
the last subsection: assign an ordered pair of subscripts 
K^Ki+i to {1 < i < M — 1), and assign KiKm to vq. 



There will then be L/2 other, special trajectories {L = 
2 for the kinetic energy, and L — A for the interaction), 
each of which starts from the diamond and ends at it. 

Case 1: L = 2, and there is one special trajectory 
vqKiVi ■ ■ ■ Kmvm (vq and vm are the two sides of the 
diamond, and Ki ~ Km)- Assign an ordered pair of 
subscripts K^Ki+i to each Vi {1 < i < M — 1). Multiply 
the whole expression by (—1). 

Case 2: L — A, and the two special trajectories are 
sKi ■ ■ ■ K\h and SK2 ■ ■ ■ K^h, where s and h are the solid 
and the hollow sides of the diamond, respectively. So 
each trajectory is the same as the one in Case 1, and 
the two (— l)'s (overall factors contributed by both tra- 
jectories) cancel. Assign an ordered pair of subscripts to 
each vertex (other than the diamond) in the same way as 
above; the diamond contributes the factor C/K^K^K'^Ki, ■ 

Case 3: L — A, and each special trajectory 
oKi • • • Kmo starts from and ends at the same side of 
the diamond. It can then be treated like an ordinary in- 
ternal trajectory; the pair KiK^f (exactly in this order) 
is filled to the first two subscript positions of [/ if o is the 
solid side, or to the last two subscript positions of [/ if o 
is the hollow side. The order of the fermion operators in 
Eqs. (121) and Q has been chosen to ensure the validity 
of this simple rule; such order has also affected the rule 
in Case 2. 



E. Thermodynamic orders and a property of the 
wave function 

The theorem derived in concerning the thermody- 
namic order (see Appendix ^ of any aforementioned 
diagram, namely 

Q = 7V3 - L/2, (9) 

is still valid, because the prerequisites of this theorem are 
unaffected by the change from a bosonic superfluid to a 
fermionic one; the fermion sign does not affect these ther- 
modynamic orders. L is the number of external points of 
the diagram. This theorem describes the scaling behavior 
of any diagram T) in the thermodynamic limit: T) ^ Vfi . 

As a corollary, any diagram in the expansion of E — fiN 
satisfies Q = 1 and is thermodynamically significant. 

In almost the same way as in ^ , we can also establish 

(v) 

a general property of the coefficients Oki - k ' namely 
it is zero whenever a nontrivial subset of fermions con- 
serve spin-momentum (that is, whenever the sum of a 
nontrivial subset of the subscripts is zero). The basic 
idea is that any thermodynamically significant contribu- 
tion to any physical observable from a coefficient whose 
nontrivial subset of subscripts conserves spin-momentum 
(called reducible coefficient) is only present at the "bot- 
tleneck" vertex of a dead end (see Appendix IB)| . but that 
dead ends can be absorbed by redefining lower order co- 
efficients. Consequently, there are no dead ends in any 
diagrams. For more details, see 
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III. FEW-BODY WAVE FUNCTIONS 

In this section we define some few-body functions and 
discuss some of their properties. They are completely 
independent of any many-body physics whatsoever. 

But the contents listed below are restricted to those 
that are directly needed in Sec. IIVI where we do a many- 
body low-density expansion. 

The internal wave function c/j-k of an isolated molecule 
satisfies {H ~ Em) J2k <^k'^k'^-kIO) = 0- 



''KK'fK' 



K' 



0, 



(10a) 

(10b) 
(10c) 



K 



where 



DkK' = (fc^ - E,n)SKK' + ^Uk,-K,K',-K' (11) 

is a hermitian matrix for the relative Dynamics of two 
fermions in opposite spin states isolated from all other 
fermions. Since Em < and the interaction is finite- 
ranged, the spatial representation of must be finite 
ranged, so 0k(T depends smoothly on k. Moreover, ob- 
viously 
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Sometimes we elect to abbreviate the spin-momentum 
Ki as i, and — Kj as i, when they do not lead to confusion, 
i = 1, 2, 3, • • • , 9. Without this notation, some equations 
would be too lengthy and cumbersome. 

The two-molecule zero-speed scattering wave 
function (^t^^T^^Ti-^T^. satisfies \H — 2E„ 



/,(4) 

/,(4) „t „t„t„t 
234 V^1234'-l''2''3'-4 



|0)=0. 



•'1234 



^2134 



'^3214 



4^4-0 (ifKi 



V^4231' 

K4 ^ 0), 



(12a) 
(12b) 



[{kf + ■■■ + ki)/2 - 2Em - e^'^ Cl>[%, + i^C/i2564634 



56 



(4) 
564 



2^^135605264 + ^14560^^36 + 2^^2356^1 

56 56 56 

+ + IT.^3456& = 0, (12c) 



56 



56 



where e^*^ ^ fl^^ is an extremely tiny positive energy, 
and is always negligible unless two of the four spin- 
momenta (Ki • • • K4) are opposite; Ref. contains a 
simpler analog. When the distance r between the centers- 
of-mass of the two molecules is large [3> max{l, r^)], the 



four-fermion wave function (in the coordinate represen- 
tation) is proportional to (1 — a,„/r) times the internal 
wave functions of two isolated molecules. Consequently, 
we have the infrared boundary condition 



47ra„ 

7V 



/'Ki^K. +0(g°) (I2d) 



when q is small or zero. Here 

, (4)q ^ , (4) 
•PKiKs — '?^Ki+q/2,-Ki+q/2,K2-q/2,-K2-q/2- 



(13) 



The global coefficient of is set by Eq. lfT2H)l . Note 

that ?!'k!k2 delta-function singularity at Ki = ±K2, 
but the intensity of the delta function approaches a finite 
constant when q ^ 0. The term 0{q^) in Eq. Ijl2dp refers 
only to the fact that the correction term approaches some 
function of K1K2 which is independent of q, when q — > 0. 
The following decomposition will be useful later: 

<^1234 = '/''1234 + Si^lhiM'i ~ '^13'52401'/'2 + (Jnf^ss^il 02 , 

(14) 

where 0K1K2K3K4 ^ ^ ^ contains no delta-function 
singularity and is also completely antisymmetric. If 
Ki + • • • + K4 = but Ki + K, 7^ (for all 2 < i < 4), 
Eq. 112c|l becomes 



[{kl + ki + kl + kl)/2 - 2Em] 0'i234 
^ ^E^125605634 + ^E ^13^6 05264 + ^E^145605 



5236 



56 



56 



56 



2E'^23560i564 + ^^^^245601536 + 0^^345601256 
56 56 56 

- £^123403 04 + C^1324 02 04 " C^1423 02 03 

- C/231401 04 + C/241301 03 " ^^341201 02 = 0. (15) 



A shorthand similar to Eq. 

0'kiK2 = 0ki+q/2,-Ki+q/2,K2-q/2,-K2-q/2- 



(16) 



Obviously ^'^'a = -(47ra,„/ilg2)0^02 + 0(g°) as g ^ 0. 
Projector onto the function space orthogonal to 0k: 



^KiK2 = '^KiK2 - 0K10K2/2- 

"Deviation functions" : 



(17) 



(18a) 



K3K4 

= rfic-^^E^KK3l0K3l'0K3, 
K3 

d'^^^l lim i ^KK3 K3K.0K., (18b) 
K3K4 

5k,k2 ^ " 1™ E ^k.K3^K2K4 0K3K4> (18c) 
K3K4 

5k,K2 ^ ^ 1™ E ^k,K3^K2K4 K3K4- (18d) 

K3K4 
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These limits exist because the divergent component in 
'/'it^*' or (j)'^2i for q ^ but q — > 0, is proportional to 
4>i4'2/q^, which is removed by the projector. d-K and 
are of the same thermodynamic order as 0k , but orthog- 
onal to 0k: Sk '/'k'^K = and X^k '^k'^k ~ 0; like- 
wise J2l <l^*l9'l2 = E2 <^2ff'l2 = El <^l5l2 = E2 ^*29l2 = 0. 

Moreover, (i_K = — ^k, and d' p, — —d'y^; g'12 = ^.9i2 ~ 

-9'i2 = 521 ; 512 = -5i2 = -5i2 = 521; g'i2 ^ No 
evidence for the vanishing of any of these functions is 
found. 

The identity 



2 



C/ll22'/'2 - -2(fc? 



E„ 



and its complex conjugate (note that Ui23i — U3412) are 
often used (explicitly or implicitly) in the next section. 
Two more identities: 



1 



A* ^(4) 



2 XI ^1234<?'2<?'i234 + 2 ^ L7l23492<?'l234 
234 234 



27ra„ 



X! ^123401 02'/'i234 = -47ra„i/fi. 



(21) 



1234 



Equation H2U|) is proved in Appendix^ its inner product 
with tpl yields Eq. (|21|l . Expressing 0^4) and d in terms 
of </)' and d', we get (to be used in Sees . II V Ul and HV H|l : 



^ [^123402 01234 + ^ X! ^123402 '^'1234 -^12|'/'2p02 
234 234 2 

+ I'/'ll^X! ^112202 - 01 X!(^1212 + Ui2l2)\(j)2\^ 
2 2 

i + hE^-4, (22) 



X! t^l234'/'l02'/'l234 + X! ^1122 I 01 P-/"! 02 
12 

- 2^(71212 101 n02 1' --47ra,„/r!. (23) 



1234 
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The three terms on the left side of Eq. (|23|l are all real. 

Using a similar method as in Appendix ^ starting 
from Eq. we can show the following identity: 

(^1234 01234 + ^1234 01234~^12340I234~^1234 01234) 

^ 34 

+ t^ll2202 + ^^221101 - (t^l212 + t^l212 + ^^1212 + t^l2l2)0102 
=^E[^13(^302+532)+^23(0l4+5l3)]- 



27ra^ 



D1(P2: 



(24) 



when Ki ^ ±K2. This will be used in Sec. UVHl 

Finally, we have the (p/2)-molecule zero-speed scatter- 
ing wave function (p = 2, 4, 6, • • • ), the generalization of 

Eqs. (nnj) and (US). 



pEn 



1=1 



2 



E 

1<2<J <p 



K.iK'Ki 



KiK,K'K" 



= (25) 



if Ki 



K„ - 0; 



/,(P) 



= otherwise. 



41. 



is antisymmetric under the interchange of any two spin- 



(19) momenta. If p = 2, e^^) =0. If p > 4, e^P) ~ 



and 

this extremely tiny energy is negligible unless the sum of 
a nontrivial subset of the p spin-momenta vanishes. 



IV. DIATOMIC MOLECULES: LOW-DENSITY 
EXPANSION 

At low density, the ground state is a dilute gas of 
molecules, any two of which are usually - but not always 
- far apart compared to the size of a molecule. In this 
section we study the ground state, exploiting such a low 
density. 



A. Low-density orders 

If the leading order contribution to a quantity x scales 
like n^n'^ in the low-density regime, we say that the 
low-density order of x is R, and may schematically write 
X ~ n^; more details are in Appendix IbI 

It is reasonable to expect that the ground state energy 
is, to the lowest order in density, simply E = EmN/2. 
So /Lt = Em/'2- + h.o.c. [h.o.c. = higher order corrections 
(Appendix El] J and the low-density order of /i is 0, in 
contrast with the case of structureless bosons 2] . 

Since /J. < 0, the low-density order of k"^ /2 — is 0, for 
any k. The low-density order of C/K1K2K3K4 is 0. 

The low-density order of Q!k is 0.5. If we only retain 
the lowest order term in the exponent in Eq. (Q, we 
get riK — |aKp/(l + I^Kp)- At low density, the spin- 
momentum distribution should approach that of fermions 
in an isolated molecule times a constant, so the width 
(~ 1 /r,„) is roughly independent of n, and the occupation 



number of each K state should be of the order 
,0.5 



< 1. 



So ckK ^ n" 

The low-density order of /3K1K2K3K4 (where Ki - • 
K4 = 0) is usually 1 (see below for exceptions), because 
|/3(RiR2R3R4)P (norm square of the Fourier transform) 
- for a given set of locations in a small region independent 
of n - is approximately associated with the probability of 
finding two molecules (or rather four fermions) in such 
a region, and this probability is approximately propor- 
tional to . 

In general, the low-density order of o^k ...k usually 
p/4. 
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The first exception to this estimate appears in 
/3K1K2K3K4) since it is approximately proportional to the 
four-fermion wave function describing the zero-velocity 
scattering of two molecules, in the low density limit. 
When the two molecules are well separated, their spa- 
tial wave function should be of the form 1 — a„i/r times 
the internal wave functions of the molecules, where r is 
the distance between the molecules. Taking the Fourier 
transformation, we see that /3K1K2K3K4 scales like amlo^ , 
when Ki -I-K2 = q ^ (the term (5q_o is encoded in ckk)- 
If there were no many-body effects, this would be an in- 
frared divergence at small q. But again, just like in the 
case of structureless bosons 0, the divergence is regu- 
lated when q ^ 1/^, where ^ is the superfluid healing 
length. We will confirm in the concrete calculations that 
^ still scales like n""'^. By reducing q = Ki + K2 from 
the order 0(1) [for any momentum k, by fc ~ 0(1) we 
mean that k has a finite and positive lower bound when 
n — > 0] to the order rfi-^ , /3K1K2K3K4 must increase from 
the low-density order to the low-density order rfi , in 
accordance with the 1/q^ dependence. After this, when 
q is further reduced, the l/q^ dependence breaks down, 
and the order of magnitude of /?KiK2K3K4 will not in- 
crease again. 

Like in the case of structureless bosons Q , we can ex- 
tend this analysis to higher orders, and arrive at a general 
result 



p/i-M+l 



(1 < M < p/2), 



(26) 



where the p spin-momenta - depending on their values - 
are divided into M l-clusters. Each 1-cluster is a collec- 
tion of (at least two) spin-momenta such that 

1. their sum is of the low-density order ^/n (if M > 2) 
or is zero (if M = 1); 

2. the sum of any nontrivial subset of them is of the 
order 0(1) [by this we mean that either the sum of 
spins does not vanish, or the sum of momenta is of 
the order 0(1)]. 

If RI > 2, the sum of spin-momenta in each 1-cluster must 
not be zero, or the whole vanishes (see Sec. Ill E)l . 

There are still exceptions to the above result, because 
of the antisymmetry of Q!^|...j^ . When the difference 

between two subscripts approaches zero, ckkI - k ^ ^■ 
However, if we regard 0(riP/"'^^^+^51^^''/^) (see Ap- 
pendix lB)) as an upper bound of a^|...j^ , the Fermi statis- 
tics does not break it. 

Now we decompose the p/2-th order (de)generator into 
the sum of a set of dispersed vertices, like in the case of a 
dilute gas of structureless bosons ^ . The decomposition 
of the first three generators are shown in Fig. |21 for the 
degenerators, we just replace each solid circle with a hol- 
low one. Despite the striking similarity between Fig. [21 
and a similar figure in 0], the corresponding coefficients 
have vastly different orders of magnitude, because of their 



-•- 

a 



\/ 
— • — 

/\ 

Y 



^0.5^ 



\/ 



\/ 



/\ 



(nO-5) +_. . ..._ (n-0-5) 
V21 ^Viii^ 



FIG. 2: The decomposition of a, /3, and 7 [the first three 
terms in the exponent in Eq. into dispersed vertices in 
the low-density regime. Each dispersed vertex is denoted by 
V with some subscript(s) indicating the number of pairs of 
fermions in each cluster; the low-density order of its associated 
coefficients in momentum space is noted. Dotted lines connect 
different clusters in v. In each cluster, the sum of spins is zero 
and the sum of momenta is of the order ^/n (if there are two 
or more clusters in v) or (if there is only one cluster), but 
for any nontrivial subset of fermions in the cluster, either the 
sum of spins does not vanish, or the sum of momenta is of the 
order 0(1). 



different thermodynamic orders. Nevertheless, their low- 
density orders are the same. Since we have proved that 
the total thermodynamic order of any diagram in the 
expansion of i? — fiN is 1 , we can concentrate on the low- 
density orders when estimating the orders of magnitude 
of such diagrams, without worrying about the thermody- 
namic orders of individual vertices. 

In such a decomposition, a = a*-^) remains a sin- 
gle term, denoted by ui. P = 

is decomposed to 

two terms, V2 and fn; the first one is restricted to 
k > kc, and the second one to A: < kc, where k — 
iiima--+aj=o,i<i<j<4:\'^i + | , and kc is a momentum 
scale satisfying 



\/nax <Cfcc < min(l//, l/r™), 
kc/\/nax independent of n. 



(27) 



where Ox is independent of n and Q; we will see that 
is actually Om- Further, 7 = a*-^-* is decomposed to three 
terms (Fig.|3), and we will see that the second term, V21, 
is most important in the low-density regime, similar to 
the situation for structureless bosons 

Such division of the spin-momentum configuration 
space is crucial for the practical computation. By decom- 
posing diagrams according to such division of the basic 
vertices, one will be able to analyze the orders of magni- 
tude of diagrams and compute their actual values much 
more easily. The physical result, however, is indepen- 
dent of such division of the spin-momentum space, and 
in particular, independent of the choice of k^ as long as 
the above inequality is satisfied. 

The low-density order of the coefficients a^^^^}^j^ ^ ^ in 
any dispersed vertex v is 

R{v) = ^-N[iv) + N;iv), (28) 
where p{v) is the number of solid lines attached to v, 
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N[{v) is the number of dotted lines contained by v, and 
N'j{v) is the number of junctions of these dotted hues. If 
M < 2 [see Eq. N[{v) = M - 1 and N'^{v) = 0; 

but if M > 3, N[{v) = M and N'Jv) = 1; see Fig. |21 
Equation (|28l) is equivalent to Eq. H26() , but will be more 
convenient in the next subsection. 



B. Low-density orders of diagrams 

To compute physical observables in the low density 
regime, it will be much more convenient to decompose 
the (de)generators in the above way first. The physical 
observables are then expanded in terms of the dispersed 
vertices, and we obtain a new set of diagrams containing 
small internal circles, solid lines, and dotted lines; for 
E — fiN we also have a diamond in each diagram, and 
for (O) {O is a product of L basic fermionic operators 
in spin- momentum space), we have L external points; in 
this paper we use the term low-density expansion in a 
narrower sense, to refer to these diagrams. 

To obtain the correct value of a physical observable 
up to a certain order in the low-density expansion, we 
must take into account all the diagrams for this observ- 
able up to this order (called significant diagrams). The 
omission of any significant diagram may invalidate the 
result completely. It is therefore crucial to derive a gen- 
eral power-counting formula, using which we can identify 
all the significant diagrams for an observable up to any 
certain order. 

Such formula is 

/ 

i? = iV3+^(P,/4-l) + 7V^/2, (29) 

i=l 

where R is the low-density order of any such diagram 
(denoted by V below), is its number of disconnected 
parts, / the number of islands. Pi the P-number of the 
i-th island, and TVj is the number of independent loops in 
the skeleton diagram. The terms are defined as follows. 

If we remove all the dotted lines (if there are any) and 
their possible junctions (see previous subsection) from 
V, but leave all the other components unchanged, 2? is 
changed to / disconnected parts, each of which is called 
an island. Different islands are either linked by dotted 
lines only, or completely disconnected. Any small circle 
(other than the external points) in an island is attached 
by a certain number of solid lines; the sum of all such 
numbers over the whole island is called the P-number of 
the island, and the sum over the whole diagram is the 
P-number of the diagram. Finally, if we reduce every 
entire island in I? to a point, but leave the dotted lines 
(if there are any) and their possible junctions unchanged 
(the dotted lines are still attached to the islands and 
junctions), we get the skeleton diagram of P. Obviously, 
if T) contains no dotted lines, its skeleton diagram is A^3 
isolated points and I = N3. 



Proof. Only two factors contribute to the low-density or- 
der R of V: 

• each dispersed vertex v contributes R{v) in accor- 
dance with Eq. if^ : 

• each independent small internal momentum q 
(^ y/n) which flows through a dotted line 
contributes 1.5, since the measure X]|q|~^ — 
O /^^^d3g/(27r)3 scales like n^-^. 

The number of independent q's above is clearly N2, as 
momentum is conserved in any part of the diagram. So 
R = P/4 - N{+ N'j + 1.5N^, where P is the P-number 
of T), N[ the total number of dotted lines, and TVj is the 
total number of junctions of the dotted lines. 

Now turn to the skeleton diagram V . It has / -I- N'^ 
vertices (islands and junctions), N[ lines (dotted hues), 
N2 independent loops, and — disconnected parts 
(the number of disconnected parts is not changed when 
the islands are reduced to points). So using Eq. IjBip . we 
get I + N!j~N[+N^-N3^ 0. 

Using the above equation to cancel the term —N[ + 
N'j in the expression of R, we get Eq. I|29|) (note that 

obviously P = ^i)- ^ 

The central result of the last subsection, Eq. H26(l . is a 
special case of Eq. H29|) . 

To better understand R, we further study Pi . Suppose 
that the contribution to Pj from all the small solid cir- 
cles (other than the external points) in the «-th island is 

Pi'^\ and that the contribution from all the small hol- 
low circles (other than the external points) in the same 

island is P/ Suppose also that the number of solid 

external points in the i-th island is L^i'^\ and that the 
number of hollow external points in the same island is 
Obviously 

i^: = + P/"\ (30a) 



P^^+> + L\+> = P^^-> + L\-> > 0, (30b) 



P}+>,P^ ^ €{0,2,4,6, 8,- ••}. (30c) 

The contribution to the balance of particle number from 
the diamond (if there is any) is cancelled. 

For convenience, we shall call any island which con- 
tains two or more external points, an external island; 
any island which contains no external point, an internal 
island. Any island which contains a diamond, a hamil- 
tonian island. Finally, any island which contains neither 
a diamond nor any external point, a (vacuum) bubble 
island, since it has the same topological structure as a 
vacuum bubble. 
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An corollary of Eq. (|30|l : for any internal island i (a 
bubble island, or a hamiltonian island which contains 
no external points), Pi can only be 4, 8, 12, 16 • • • , and 
Pi/A — 1 can only be a nonnegative integer. It is this 
property that makes Eq. H29|l particularly useful in the 
search of all the diagrams in the low-density expansion 
of an observable up to a certain order. 

The bubble island whose P-number is 4 has a unique 
topological struture - a solid circle linked to a hollow one 
by two solid lines - and contributes nothing to R; it will 
be called the simplest bubble. 

An alternative (sometimes useful) form of Eq. 129|) is 

i? = iVg + (-Pext/4 - /cxt) + + (Pint/4 - /int), (31) 

where I^xt and /int are the numbers of external islands 
and internal ones, respectively; Poxt and Pint are the total 
contributions to P from the external islands and internal 
ones, respectively. Pint /4 — /int is always a nonnegative 
integer. 

In summary, V is more conveniently viewed as an 
"archipelago" (islands with weak links) when its low- 
density order is analyzed. The archipelago structure is 
a direct consequence of the presence of two sets of well- 
separated length scales in the problem: the superfluid 
healing length (x l/y/n, and the few-body length scales I, 

; etc. 

When computing its actual value, we of course view V 
in the "usual" way, namely as a set of dispersed vertices 
with a diamond or/and external points. 

C. E — fiN to the order 

In any diagram in the low-density expansion of E—fiN, 
there can not be any "dangling" dotted line (one whose 
removal would change such a diagram to two discon- 
nected parts), since dead ends are absent fSec. lIIEll . So 
if there is any dotted line in such a diagram, iVj > 1 and 
R > 1.5 [see Eq. (P|)]. 

So if P = 1, there can only be one island, the hamilto- 
nian island. Moreover, the P-number of this island must 
be 4. So we can only find two diagrams satisfying P = 1, 
in the expansion of P — /iiV. They are Ti and T2 shown 
in Fig. |S1 

^1 = II(^V2-m)I«k|', (32a) 

K 

^2 = ^ X! ^K,-K,K' -K'ttKaK'- (32b) 
KK' 

We then adjust the parameters to minimize E — fj.N: 
d{E - fiN)/da*j^ = 0. So 

fc^^K + ^ X! ^K,-K,K',-K'"K' = 2/iaK + h.o.c, (33) 
K' 

where h.o.c. stands for higher order corrections. Since the 
fermions have formed bound pairs, and the density is low. 



ttK should be approximately proportional to the internal 
wave function of an isolated molecule, and 2fj, « Em- 
Comparing the above equation with Eq. 110a|l . we find 
Q!K = ri4>K + h.o.c. The number of fermions is iV = 
EkI"kP + /1.0.C. = 2\r]\'^ + h.o.c, so 77 = y^N/2 + h.o.c, 
where we have chosen the phase factor of 77. We thus get 



ttK = \/N^(l)K + h.o.c, (34) 

fi:^ Em /2 + h.o.c, (35) 

E/n:^nE„j2 + h.o.c, (36) 

where E is obtained by solving dE/dN — /i, and 

N^ = N^=N/2, (37a) 

n| = n| = n/2 = N/{2n). (37b) 



D. First cluster-separation theorem 

In this subsection and the next, we discuss an elemen- 
tary property of the system in the low-density regime. It 
will help to further the low-density expansion. 

Consider p spin-momenta Ki • • • Kp, where p is posi- 
tive and even, Ki -I- • • • -I- Kp — 0, and the sum of any 
nontrivial subset of these spin-momenta is nonzero. Sup- 
pose that they are divided into M 1-clusters (M > 1), 
and the I'-th (1 < < M) 1-cluster contains Pi, spin- 
momenta: Kfc,^+i, • • • ,Ke„, where YltLiPi^ = Pi 
61/ = P1+P2 + - ■ ■+Pv, and bi, — e^~p^. Suppose also that 
there does not exist any 1-cluster other than these M. 
Make the decomposition Tib„+i — '^i,u + <\v/Pv (1 < « < 
Pu), such that Y^^i=i "^i^v = 0- Consider the ground state 
expectation value F^^ i> = (c~ c~ ■ ■ ■ c~ ). The 

Ki---Kp \ Kp Kp_i Ki' 

first cluster-separation theorem states that there exists 
some coefficient rj, which depends on qi • • • but not on 
KijKa,! • • • Kpi_i, • • • , Ki3,/K2,A/ • • • Kp^^^^M, such that 

M 

j^k,...kp-^n<iK.,„...Kp„,„ (38) 

with a relative error ~ . 

There are infinitely many diagrams in the low-density 
expansion of F , but their low-density orders have a lower 
bound. Consider any one of these diagrams, V: N3 > 1, 
Pcxt > P, and /cxt < M (it is impossible to separate the 
p external points into more than M islands, since spin- 
momentum is conserved at every internal vertex, and the 
dotted lines can not carry spin or any large momentum) , 
so we deduce from Eq. that R{V) > p/i - M +1. 
The low-density order of F is therefore P = p/4 — A/ + 1 
or higher. One of the leading contributions (with low- 
density order P) to P is ~ . 

While it is intuitively clear that the right side of 
Eq. (|^ should be the form of the lowest order approx- 
imation of P, it is not a priori obvious why the next- 
to-leading order correction (of the order ti^+O-^) jg qI^q 
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proportional to such a product of zero-speed scattering 
wave functions. 

The proof is quite involved. In Appendix ICl we prove 
this theorem for the cases in which each 1-cluster contains 
only two spin-momenta. To compute E — fiN up to the 
order n^-^, such an incomplete proof will be sufficient. 
The complete proof will appear in a future work. 



E. Second cluster-separation theorem 

There is a similar formula for the coefhcient of the 
dispersed vertex: 

M 

4:...K.='^'n<iK,....K„.. (39) 
1^=1 

plus a correction term whose low-density order is i? + 1 
or higher, where R = p/4— M + 1, and the q's, K's, and 
K's are defined in Sec. lIVDl The coefhcient rj' depends 
on qiq2 • • • (Im and n, but is independent of the K's. 

The nontrivial content of this theorem is that the next- 
to-leading order correction (~ ^^le coefficient 
of the dispersed vertex is factorizable in the same manner 
as the leading term. 

In Appendix ^ we give the proof of this theorem for 
the cases in which each 1-cluster contains only two spin- 
momenta; it builds upon the results already proved in 
Appendix [HI and is even more involved. For the deter- 
mination of E — up to the order n^-^, such special 
results will suffice. The complete proof will appear later. 

This partial proof will simplify the logic in the concrete 
calculations of E — fj,N. More importantly, it enables us 
to accurately determine some other interesting quanti- 
ties, such as the many-body effects in the momentum 
distribution of the fermions. 

The second cluster-separation theorem is more fun- 
damental in the present formalism, since the first one 
and some other similar properties, such as that of 
(4i+q/2^-Ki+q/2^-K2+q/2CK.+q/2). ^re its corollarics. 
The first one is presented for two reasons: 1) it con- 
cerns physical observables rather than parameters in a 
particular calculational framework, and 2) the proof of 
the second theorem is built upon the proof of the first. 
Hereafter we use the name "cluster-separation theorem" 
to refer to the second theorem. 



F. Diagrams for E — fiN up to the order n^''' 

In this subsection we identify all the diagrams in the 
low-density expansion of E — fiN up to the order n^-^. 

Any such diagram V has — 1 since it is a connected 
one. So i? < 2.5 implies < S [Eq. 

If N2 = 0, there is no dotted line at all in V. So there is 
only one island, the hamiltonian island, whose P-number 
can only be 4 (i? = 1) or 8 (i? = 2). This is because 




FIG. 3: A pair of diagrams Di and 1)2', their sum is smaller 
than each one of them by a factor of the order . Their 
structures within the "X" -circles are identical, namely a cer- 
tain number of bubble islands of certain structures, linked by 
dotted lines in a certain way. 

any island in the expansion of E ~ iiN is an internal one, 
whose P-number must be a multiple of 4 fSec lIVB)) . The 
two R = 1 diagrams are Ti and T2 (Fig.jSJ, whose sum 
however has a higher low-density order, R — 2. This is 
the first instance of the "R+V^ -cancellation mechanism, 
discussed in the following. 

Consider an arbitrary pair of diagrams in the low- 
density expansion of E — /.iN, Vi and I?2, whose hamilto- 
nian islands are topologically the same as Ti and T2, re- 
spectively: if their other parts are identical, we call them 
an "i? + 1" -pair. The two diagrams obviously have the 
same low-density order, R. The second cluster-separation 
theorem implies that the low-density order of their sum 
is at least R-\- 1. 

If Vi and I?2 are just Ti and T2, the proof is sim- 
ple: substituting an = + 0{n^-^) (Sec. IIV E|l into 
Eq. (I22Jl, and noting that /i = E,n/2 + 0{n^) (Ap- 
pendix Q, we immediately get this result. 

If Vi and I?2 contain dotted lines, they can be written 
in the form [see Fig. 01 and Eq. (|77|) ] 

^1= II (fcV2-S„,/2-/I + gV8)^^K. (40a) 

9<fee;K 

9<fec;KiK2 

where Jl = fj, — £'m/2, the first subscript of X^^j^^ is 
associated with the creation of a pair of fermions with 
spin-momenta K2 + q/2 and — K2 + q/2, and the second 
subscript is associated with the annihilation of a pair 
with spin-momenta Ki 4- q/2 and — Ki + q/2. In Vi, 
an average over K and — K has been taken, since X^^j^^ 
is antisymmetric for K2 — K2 (and similarly for Ki), 
and thus is symmetric for K — > — K. In 7^2, U has 

Galilean symmetry and is independent of q fSec. Ill All . 
Note finally that the signs of Vi and 2?2 are correct, as 
is clear if one restrict every trajectory (Sees. Ill CI and 
IIID(I within an island when determining the sign. If the 
small solid circle associated with the first subscript is 
contained by a dispersed vertex • • • 1 ' ^^^^ from our 

A/ 

established result we find that X^^^^ — (f>K2X'^_^ with 
relative error of the order n^, and the R -\- 1 cancel- 
lation between Vi and I?2 immediately follows. If the 
small hollow circle associated with the second subscript 
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FIG. 4: Two diagrams with R — 2.5 each. Their sum should 
be ~ n"^'^, but the proof is not yet given. However, trivially 
this sum is at least of the low-density order n^. 

is contained by v^^ _ then ^KaKi — '?^ki^"K2 ''^ith 

M 

relative error of the order n^, and the R + 1 cancella- 
tion follows from the complex conjugate of Eq. (|l()a(l . 
For R < 2.5, there is only one "R + l"-pair for which 
we have not shown the R + 1 cancellation, purely due 
to our incomplete proof of the cluster-separation theo- 
rem. This case is shown in Fig. 0] It involves the vertex 
V21 and V21, and each diagram has R — 2.5. However, 

/3Ki+q/4,K2-fq/4,K3+q/4,K4+q/4,K'-q/2,-K'-q/2 (Q V^: 

Ki -I- • • • -I- K4 = 0) is proportional to (/)k' , at the leading 
order at least, since it is approximately associated with 
the correlation of four fermions and a molecule which is 
far away (with distance ~ ^) and the internal structure 
of that molecule is almost not influenced by the first four 
fermions. So the sum of these two diagrams is at least 
of the low-density order n^, and can be omitted in the 
computation of — fj,N up to the order n^-^. 

The counterpart of the "R-l-1" mechanism in the con- 
text of structureless bosons |3| is a much more obvious 
fact that — /Xboson is of the order when the boson 
momentum g = or g ~ n^/^. 

Because of the "i?-|-l"-cancellation mechanism, the set 
of significant diagrams, up to the order n^-^, is greatly 
reduced. If the "i? -I- l"-pair Vi and I?2 satisfy R > 
2 each, then their sum is at least of the order and 
negligible. 

If the hamiltonian island is Ti or T2, only the N2 — 
or A'^ = 1 diagrams are significant, and the bubble 
islands are all the simplest ones (because of the R + 1 
mechanism). They are Ti, T2, Tia and T2a, shown in 
Fig.El 

If the hamiltonian island has P = 8, then N2 = {R — 
2) or A^^ = 1 (i? = 2.5). If A^^^ = 0, the hamiltonian island 
(the whole diagram) may contain two Vi^s and two uj's 
(Ta, T4, T5, and Tg), or one V2 and two uj's (Ty), or two 
vi's and one (^s); or one V2 and one (Tg and Tio). 
li N2 — 1, we just exhaust all the topologically distinct 
ways of attaching a dashed line to the hamiltonian island, 
and obtain the diagrams T^a, Tafc, ^Sc, Tka, 74b, T^c, T^d, 

TI rf\ rri rji rr~i rri rr\ rri rT~i rri rri rri 

5a, J5b, -^5c, -^Sd, J^6a, J^6b, ^60, J6d, ^Ta, ^76, ^Sa, ^86, 

Tga, and Tioa- The bubble islands in these diagrams are 
all the simplest ones or i? > 3.5 [Eq. (HHH- AH these 
diagrams are shown in Fig. [S] 

Each dashed line represents a geometric series associ- 
ated with a variable number of simplest bubbles. Each 
end of the dashed line is attached by an arbitrary posi- 
tive even number of solid lines. The concrete meanings 
of the dashed lines used in Fig.|Slare shown in Fig. El 
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FIG. 5: Diagrams in the low-density expansion of i5 — fiN 
of molecules (pairs of fermions) up to the order n^'^ . The 
composite vertices (each of which contains a dashed line) are 
defined in Fig. |H] each of them represents a geometric series. 
In Tsb, Tzc, Tid, and T^d, a cross is marked on the dashed line, 
to indicate the absence of the first term in such a series, ie. 
/3K1K2K3K4 (or its complex conjugate), because two of these 
four spin-momenta happen to be exactly opposite. 

In Tab, Tsc, T^d, and T^d, a cross is marked on the 
dashed line, to indicate the absence of the first term in the 
geometric series, ie. /3K1K2K3K4 or its complex conjugate 
(recall that a diagram can not contain any dead end). 

G. E - fiN to the order n'^ 

The diagrams for E — fiN up to the order are Ti 
through Tio only. Note that Tia and T2a are an "R + 1"- 
pair and their sum is of the low-density order n^'^. 

Ti =^(fcV2-Ai)|«K|', (41a) 

K 
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FIG. 6: Definition of some composite vertices, each of which is a geometric series. The spin-momenta K's are of the order 1, 
and the momentum q ~ y^. The dispersed vertices are defined in Fig. |21 



^2 = 4 H t^Ki,-Ki,K2,-K2"Ki"K2' (41b) 

K1K2 

^3 = -E(^V2-i^™/2)|aK|^ (41c) 

K 

^4 = ^ C^Ki -Ki,K2.-K2aKil"K2paK2, (41d) 

K1K2 

^5 = -^ X! ^Ki -K1.K2 -K2|aKiP"ki"K2' (41e) 
K1K2 

?6 = ^ X! f^KiK2KiK2|aKiHaK2l^: (41f) 
K1K2 

1 v^' 

T7 = -4/^ ?^-Ki -K2,K3K4aKi"K2/3KiK2K3K4, 
K's 

(41g) 

1 v^' 

^8 = ^4^^ t^KiK2.-K3,-K4/3KiK2K3K4'^K3Q^K4' 
K's 

(41h) 

^9 = ^E'(^?/2-^™/2)I/3kiK2K3K4|', (41i) 

K's 

TlO = t^KlK2K5K6/3KiK2K3K4/5K5K6K3K4I (41j) 

K's 

where excludes cases in which the sum of any two 
subscripts of /3 (or /?*) is less than fee, and in T3 (and Tg) 
the term containing /I is omitted since it is ~ n"^. 

The parameters a and (3 are then adjusted to minimize 
E — fiN. We first write the terms containing f3* {ie Tg + 



Tg +T10) in the form Ek's /5kiK2K3K4^KiK2K3K4; only 
the completely antisymmetric component of 9 contributes 
to the sum, and this component must vanish since it is 
independent of /3* and d{E - iJ,N)/dp* = 0. We get 

[{kf + ki + kl + kl)/2 - 2£;„]/3i234 

+ C^1256P5634 + 1^2^ ^^1356/35264 + ^^1456/35236 

56 56 56 

1^' Iv^' 1^^' 

+ 2^-^ C^2356/3l564 + 2Z_^ 1^2456/^1536 + 2/-^ C^3456/3l256 
56 56 56 

- t^l234a3a4 + C^l324a2a4 - C/l423"2a3 

- C^23i4"ia4 + C/24i3"ia3 - U^ii2aia2 = 0, (42) 

The excluded regions in the spin-momentum space for 
the summation X)56 contribute terms of the low-density 
order n^'^, which are smaller than the dominant terms by 
a factor of the order 'nP^^\ a-KS lowest order expression is 
given by Eq. So the comparison between Eq. 

and Eq. ((T5|) yields 

/3k,K2K3K4 = ^T</'k,K2K3K4 + 0{n^-^^'^) (43) 

if mincr,+crj=o,i<j<i<4|Kj -f- Kj| > k^. 

Now turn to d{E - fiN) /da*^ = 0. Write 5{E- fiN) = 
X^K ^k'^'^k + ■ ■ ■ • Only the antisymmetric component 
of contributes to the sum, since a_K = — Q^k- This 



14 



component must vanish, and we get (after a little algebra) 

-01202 — C?\ ^ D2\Oi*2 + ^ £'i2|a2|^a2 
2 2 2 

+ |aiPX!^ii22"2 -ai^(C/i2i2 + C^i2i2)|a2p 

2 2 

+ C^l234a2/^1234 + C^1234a;2/3l234 (44) 

234 234 

plus correction terms ~ n^, where Dn = Z3i2 — 2/i(5i2 and 
Jl = jj, — Em/'2. The terms on the right side of Eq. H44(l 
are ~ n^'^, and in this subsection they can be simplified 
with the lowest order formulas [Eqs. H34() and H43() ]: we 
then solve the equation to obtain more accurate results 
of a-K and /i. 

The first term on the right side of Eq. H44|l vanishes 
at the order n^'^, and the remaining terms are further 
simplified with Eq. H22|) (the excluded regions in the spin- 
momentum space for X]234 contribute terms of the 
low-density order and are unimportant here): 

^3/2 

2 2 

(45) 

We decompose aK into a component parallel to (j)-K and 
one orthogonal to 

CtK = Va't'K + ctk^ X!'^k"k=07 (46) 
K 

where ?7q is elected as positive, since there is a gauge sym- 
metry: when a^P) in Eq. Q is changed by a factor e'''^ 
{9 is real and independent of p) , the particle number and 
the energy expectation values are invariant. Substituting 
Eq. (gni into Eq. and noting that rja = y/N^+h.o.c. 
('Sec UvTl . we get 

Di2{a^ - N^^^dyn) = 2]lai + 0{n^). (47) 

2 

So 

A* = £;™/2-K7ra„n|-HO(ni-^), (48) 

and the first term on the right side of Eq. H47|l is ~ n^-^ 
and negligible; noting also that the spectrum of D12 
(when restricted to the subspace orthogonal to (/)k) has 
a lower bound \Em\, we get 

= n]''^Vt^/'^d'^ + 0(n2). (49) 

Solving di?/dA^ = /i, we get 

£;/r2 = +7ra„n^ + 0(^2-5). (50) 

As a consistency check, we may use the above results for 
a and (3 to calculate Ti + • • • + Tio, and verify that it is 
equal to -E — fiN — ~Tra„in^Q, + h.o.c. 



Firstly, Ti +T2 = -mEkI"k|' + 5 E12 al^i2a2, but 
X;i2ai^i2a2 = E12 "i"*-Di2a^ ~ n^, so Ti + T2 
— 27ram'^l^^ + h.o.c. 

Secondly, T3 + ■ ■ ■ + Tt = Nf[t(^'^ + i^^)] + 0{n^-^n), 
where - -(1/2) 0*i?i2|02|'</'2 = 0, and t(2) = 

-(l/4)El2C^1122l'/'lP0I02 + (l/2)El2C^1212|</'in02p- 

(1/4) J2 C^i2340I02</''i234 = +7ra„/^^ [see Eq. so 

1234 

+ ■ ■ ■ + Tf = +7ra„jn|f2 + h.o.c. 

Thirdly, Ts + Tg + Tw = 

(7V2/24) E,234 0':234[-6C^12340304 + 2(fc2 - E„,)^[^^^ + 

3 S56 ^i256</'5634] + ©(n^'^il), but thc antisymmetric 
component of the expression in the brackets is just the 
left side of Eq. so Tg + Jg -I- Tw = 0{n^-^n). 

So Ti + • • • + Tio = —iramn^^ + h.o.c. 

Equation (|48|l is indeed the well-known equation of 
state at the mean-field level: 

E„i + 47ra„in„/TO,„ -|- h.o.c, 

where pm — 2/i is the molecular chemical potential, = 
n| is the molecular density, and mm = 2m = 2 is the 
molecular mass. We have thus successfully constructed a 
many-body theory which is completely compatible with 
the exact, nonperturbative solution to the quantum four- 
fermion problem, with fairly arbitrary finite-range inter- 
particle interaction [not restricted to a contact interac- 
tion for which / ^ 0). 

Although common wisdom has had some success on the 
equation of state, a working many-body theory is needed 
for the understanding of many other important features 
of the system. Other predictions of this theory will be 
presented in Sec.0 Here we just point out that Eq. H49|l 
is directly related to the deviation of (c_kj,Ck|) from the 
internal wave function of an isolated molecule. For a two- 
component Fermi gas near a Feshbach resonance, such 
deviation is most dramatic in the unitarity limit, in which 
the width of the internal wave function of an isolated 
molecule in momentum space l/r„i 0, but (c_kxCk|) 
remains a finite width ^ kp ^ n^l'^ . The present theory 
helps to bridge the two regimes, and describes accurately 
how the physical picture of the system begins to change, 
as we start from the BEG limit, going toward thc unitary 
regime. 

In thc next subsection, we go beyond mean-field. 



H. E - p.N to the order ri}-^ 

E — ^iN is equal to Ti + ■ ■ • + Tio plus the following 
T's, up to the order v?^^ . 

Tia + T2a^\Y. [^12 + [q^ I ^ " 2/1)^12] C^i, (51a) 

12, q 

T^a = -2^(fci - Em)\a2?C^2. (51b) 

2,q 
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Tsb + T^c — — 



+ c.c. 



^12, q 

T4h + Tsb = -i ?^1122l"2pC^i + C.C. 

^12, q 

74c + = - ^^1122a2^j2 + C-C. y 

Tid + T^d = [ ^ t^ii22aia2^2; 

^12, q 

Tea = ^1212 I ai I 
12, q 

Teb = ^ C/i2i2Q:ia2Ci'*2 

12, q 

Tqc + Ted = 2 X! t^i2i2aia2^?2 + 

12, q 

T7a+Ts,a = - X! ^^1234^^2/^1234 + C.C 

^1234, q 

1 

T-b = 2^ t^l234a27l234;5^15> 

12345, q 

= "7 X! ^1234a37m4;5^45' 



iq 

^22' 



12' 



C.C. 



(51c) 

(51d) 

(51e) 

(51f) 

(51g) 
(51h) 
(51i) 

(51j) 
(51k) 
(511) 
(51ni) 



12345, q 



Tga = ^ XI (^l/2 ~ ^™/2)7l*234;77l234;i 



123478, q 



X [(l-/3t/3q/4)^i]^^, (51n) 

TlOa = ^ E'f^l2567l2^34;775V8[(l " ^l^^l ^TXv 

(51o) 



12345678, q 



where the arabic-number sw&scripts (except those of T) 
are spin-momenta, the summation over q (7^ 0) is re- 
stricted to q < fee, corrections of order rv^ or higher are 
omitted, c.c. = complex conjugate, 

q _ 
7KiK2K3K4;K5 = 7Ki+q/4,-- ,K4+q/4,K5-q/2,-K5-q/2> 

(52) 

excludes cases in which the sum of any two spin- 
momenta with opposite spins in the same 4-fermion l- 
cluster might be smaller than fee, and the matrices 

C7q= (l-/3q/3t/4)"'/3q/3V2> (53a) 

= (1 - /3q/3t/4)"'/3q = /3q ( 1 - /3t/3q/4)"\ (53b) 
Iq = (1 - /3q/3t/4)-'/3q/3t/3q/4 = Aq - /?q. (53c) 

For any matrix Afq, (Mq)KiK2 = ^^KiKj- Equa- 
tion ljDl|) defines /3q. 



Like in Ref. we first determine 7i'234.5, then /3^2j 
and then /3i234, and finally determine ai. 

We first express d{E — fiN) / d-fl234-7 = to the leading 
order, using the lowest order formulas q:k = \/N^<Pk + 
h.o.c, and 

/5i2 = ^q<?^i<^2 + /i-o.c, Xq = a;_q 

in accordance with the cluster-separation theorem, where 
Xq is so-far an unknown number; the second equality 
follows from the symmetry f]^2 = P-zi^ [Eq. l|Dl|l ]. So 



^12 



E3/3?3[(l-^q)'']32=a;q0iE3 03[(l-rq)-i] 



h.o.c. [vq is defined by Eq. (|D5|) ]. and 



32 - 



^ { - (Xq v/]Vf /4) t/i2340304<^8 + (1/24) (/c2-£;„)7^234;8 
8 

+ (l/16)5]Vi25675^634.8},„j(l -^q)"']87 = 0, 
56 

where {• • • }i234 means antisymmetrization with respect 
to the subscripts 1234. Since (1 — r^)^^ is invertible, 

|-6(2a;qy]Vf07)L/i2340304 + 2{kl - £'m)7l'234;7 

+ 35]Vl25675V;7}^23^ = 0. (54) 



Comparing this equation with Eq. (|15|l , we find 

7l234:5 = 2y]V|Xq05(/)i234 -|- h.O.C. 



(55) 



Next we determine (3^2 (0 < 1 < ^c)- We com- 
pute d{E — ^N)/dPl2 to the order , and then 
determine what value /3^2 should take, in order for this 
partial derivative to vanish. For this purpose ofK (in 
Tsa, ■ ■ ■ ,Tioa), 7i234;5 (i^ T76---Tioa), and /3i234 (in 
Tja + Tga) can be retained to the lowest order. After 
taking the partial derivative, we can retain /3^2 and PI2 
in d{T3a + • • • + Tioa)/dPl2 to the lowest order, since 



this partial derivative is 



Note, however, that 



d{Tia + T2a)/dPi2 IS ~ 'nP and should be treated exactly 
at this point. Note also that Tg, Tg, Tio do not contribute 
to d[E — ijlN) / d I3l2 , since in r8,T9,Tio, the parameters 
/^i234 refer to a different region of the spin-momentum 
configuration space than /3^2 (? < ^c); similarly the fac- 
tor /3i234 in is independent of the partial derivative. 

9(T3a/2 + Tia)/dPl^ = up to the order , since 
El (t^lD^2 = 0; similarly 9(r3a/2 + T^aJ/dPl^ = 0, and 
d[{T3b + Tse) + {T4d + T5d)]/d(3l^ = 0. Using the above 
result of 7, writing Aq = Pq{l — rq)~^ and approximating 
the left factor /3q with its lowest order expression, we 
easily see that 9(Tgf, -I- Tga -I- Tioa)/(?/?i2 = up to the 
order . For the remaining terms, we use the identities 



dAi/dp;^'^Xq ^ dAi^/dp;^ ^xg, 

dCq/dP;f = (l/2)/3qXq, 
dAq/dP;^' = (l/4)/?qXq/3q, 
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where the dependence of Xq on KgKgq' is suppressed 
for brevity. We get 

d{E - f,N)/dp*J = S^,X^, + 0{n'-'n-'), 

56, q 



£)?2 = (fc? - Em + q'l^ - 2-11)512 + (l/2)C/il22, 



^bQ = ~4 E ^5612 '/'5612 + 2 
12 



5234 V^2 V^5234: 



234 



^^"^ = "1^ E f^ll22l01P0I02 + J 5] C/1212|0102P 



^ E ^3412'/'1' 



12 

7ra^ 



'2<P 1234 



1234 



417 ' 



where S'^^-' is simphfied with Eq. (|23|l in the end; the sum 
of any two spin-momenta in the function (p' is restricted 
to nonzero. From the definition of ylq, we can easily 
show that EE -X^^ = -X^^ = X^^^; in S'^g, only the 
component with this same symmetry contributes to the 
partial derivative, and this component must vanish, since 

= (1 - p%/4)-\df3i/d/3;^'){i - f3qPi/4r\ 

where the matrices (1 - /?^/?q/4)^^ and (1 - /3q/3^/4)^^ 
are invertible. The resultant equation is then simplified 
with the identities in Sec. Illll in particular Eqs. (|22|l . 
(Uni, and The result is 

E ^5iPw + E ^61 /^5i + 2™™nt (1 + Axq + xD^r^^e 
1 1 

- n|(l + 2a:q) [{Dd')rM + MDd')e] 

- "T E ^515^6 -"tE ^61551 = 0, (56) 



with an absolute error of the low-density order n}^^. Here 
{Dd')5 = J2i D5id[. We then do an orthogonal decom- 
position: 

= XqM2 + Pi'^^'h + 01 /3f^^ + /3^2^^ (57a) 
^l(3['^^ ^ Y 0^/3^2^'^ ^ E ^ 0, (57b) 



where q = — q. The third term on the right side of 
Eq. 1)5 7a|) is related to the second by the symmetry 
I3i2 = ^21- Substituting Eq. (|^7a| into Eq. we 
get four equations, for the four mutually orthogonal sub- 
spaces. The linear kernal in Eq. (|56|l (the sum of the 
two D'^^s) has a single eigenvalue of the order (ie 



— 4^), associated with the eigenvector (j)i(j)2, and 
all the other eigenvalues are at least about \Em\- In the 
subspaces orthogonal to 4>i4>2, the difference between D"^ 
and D is negligible. In the subspace containing (f)i(j)2, we 
may approximate Jl with its lowest order value, Trom'^t 
[Eq. (|48() ]. since at the current order, Xq can only be 
determined to the leading order. We thus get 



(g^/2 + 47ra„in|)a;q -|- 27ra„in|(l 4- x^) = 0, 
/^Wq = n|(l -I- 2a;q)dK + h.o.c, 
^k!k2 ^n^9K,K2+h.o.c.. 
The solution to the first equation is 



(58a) 
(58b) 
(58c) 



= - (l + fq^ - eg V2 + eV) , e = (8^a„,nt)-i/2. 

(59) 

The solution with norm greater than 1 is unphysical 
and has been discarded, like in the case of structureless 
bosons (cf 0, 13) • The other components of (3i2 remind 
us that additional features are present, in contrast with 
a gas of structureless bosons. 

Now we turn to P1234 [excluding /?^2 (? < ^c)]- Its 
leading order expression is determined in Sec. IIV Gl and 
we now proceed to the next order (~ n^'^). Only one 
additional term contributes to d{E — f^N)/df3l234 — Oj 
namely Tsa- 



2{kl - i?m)/3i234 + 3^ [^1256/35634 - 6C/i234a3a4 
56 

-6C/1234 E ^34) =0, (60) 

0<9<fee ^ 1234 



where excludes the region occupied by /3q {q < kc), 
so the last term is in some sense complementary to the 
second. Since our accuracy goal here is n^-^, A'^^ can be 
retained to the lowest order, (j)i(j)2Xti/ (l — Sq), and an 
should be retained to the next-to-leading order, q;k = 
VaipK (the component orthogonal to 0k is ^ n^'^ and is 
negligible here, since it is multiplied by the other q;k' ^ 
nO-5). 

The lowest order formula /3i234 = -^T'?^i234 
combined with Eq. (|12d|) (the sum of any 
two spin-momenta is nonzero here) leads to 
/3-K3+q,-K4-q,K3,K4 ~ -(47ramn|/g2)(^304 and 
/3-K4+q,-K3-q,K3,K4 ~ (47ra„n|/g2)<^304 for 
^a,„ri| ^ g ^ min(l/Z, l/r™). If the term con- 
taining followed this asymptotic formula for all 
g < fee, it would not alter the solution /3i234 = 'ila't''i234 
at the order n^'^. It is the nonzero difference, namely 



34 - 0304(-47ra™n|/(j^) = <j)3(j)4 2^^) ' 



that constitutes an additional driving term (in addition 
to the term containing a's). Since this term has the same 
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structure as the a term, it can be absorbed into the a 
term, thus yielding the solution 



51234 



E T 



0<g</£c 



1 



^1234; 



accurate up to the low-density order n^'^. Substituting 
Eq. H59() . and noting that kc 3> y/amnf, we get 



1234 



^1234- 



(61) 



Now turn to d{E — iiN)/dal^ = 0. This equation is 
determined to the order n^-^ in Eq. (|44|l . Here we proceed 
to the order by including contributions from Eq. ()51|l . 

2 

0<q<ka ^ 

+ E ^1122 [«2(C2^2 + C^l) + «lC21 + «t A^i +«;I^2/2] 
2 

+ E t^221l(aiCi^2 + «;^n/2) + E'^l2347?234;5^;5/2 
2 2345 

- ^1212 [2aiC2^2 + 2a2Ci^2 + a2(^?2 + ^21)] j , 
9 J 11 



where {■ • • }ii retains only the component antisymmetric 
under 1 <-^- 1, and f{ ' are the terms on the right side 
of Eq. H44(l . Our accuracy goal for this equation is n? . 
Each term in f^'^'' is ~ n^'^ and is computed to the next- 
to-leading order, at which ok = ??a</'K (the components 
orthogonal to (/)k are negligible here) and /3i234 is given 
by Eq. (|61|l . Each term in z}^-* is ^ and only computed 
to the leading order, at which ckk = \/ 0k , C12 — 

(1^1^2X1 /{l-X^), = 0102X,/(1-X2), A^, ^ '\>\x\l(\- 

x^), and 7i234:5 ^ 2y]V|a;,(/)50'i234 [Eq. ((23)]. Here is 
given by Eq. . 



/I'-" = Ui[-(fc?-i?m)l01P01+E'^ll22l'^2p</'2/2 

2 

- 201 ^ C/l212|02|'l + {^l + 3i)^Vi23402 0'l234 



234 



11 

(62) 



accurate up to the low-density order n^, where 

and excludes three regions in the spin-momentum 
configuration space: IK2 -I- Ki| < /cc, IK3 -t- Ki| < fee, 
or IK4 -I- Ki| < kc (the intersection of these regions is 



smaller by a factor cx k\ and negligible), ie, 

E ^123402 01234 = E ^12340201234 
234 234 

- E f^-Ki,Ki-q,K5-q/2,-K5-q/2 
5;|q|<fee 

^ 0-Ki-|-q0Ki,-Ki+q,K5-q/2,-K5-q/2 

— (the other two small regions). 

Since 0'i234 ^ l/?^ at small q fSec. the contribu- 
tions of these regions to /}^'^'' scales like n^-^kc and is 
~ n^, if we let kc be ^/na^ times a fixed dimension- 
less coefficient (^ 1). So we only need to compute 
them at the leading order; for example, the second term 
on the right side of the above equation is simplified as 
+ 101 P E5 C^ii5505 Eo<9<fc. 47ra™/ng2. 

We can thus change Y^' to ^ in the last term of 
Eq. H62|) . and in the same time add the following cor- 
rection term (~ n^) to the right side of Eq. (|62|l : 



^-iPE ^112202 - 201 ^t/1212 102 n E 02 

2 2 ''llq<fce ^ 

(63) 



which has been retained to the leading order. On the 
other hand, the terms containing in are also 
dependent of fee, since \/^X]<;<fc^ ^12 = 0i02(3t — 
\/^'}2,q<k '^'^'^mTi^ / Q^) ■ Direct calculation immedi- 
ately reveals that the sum of the terms cx J2q<k V?^ 
in is opposite to Eq. /[^-'"^ + is thus in- 

dependent of kc at the order (the in f^^"^ can be 
simplified as ^ since their difference is of the order n^-^). 
This independence is not a coincidence, since different 
terms in this pcrturbative expansion refer to complemen- 
tary parts (in spin-momentum configuration space) of an 
underlying seamless whole. The division between q < kc 
and q > kc is just for the convenience of concrete calcula- 
tions, and the physical outcome is of course independent 
of kc- 

After straightforward calculation, we get 



fr'+fi''-{vi+^t)s- 



I6t 



^^i^l2|02p02 + O(n2-5), 



where s is exactly equal to the left side of Eq. H22|l . Thus 

- 2Jlria(j>i + ^ i5i2a^ = -01 ('7^ + 5i)27ra™/r2 
2 

+ E [ivl + 5t)dyn - (16V5)|02|'02] + 0{n^-^). 



The difference between ^2 ^12(^2 ''^^'^ Y2 ^12.(^2 
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n and negligible. Noting r]o 



/]V| + h.o.c, we get 



-K K 



^^K ^ 



l + (40/3V^)(n|a^)^/^ tto m/^ + 0{n ) 



-K K -K K -K 

o • o • o o • o • o o • 



K 



(64a) 



l + iA0/3V^){n^alf^]d[/n 
- (128/15V^)(nTa^)'/'iV^'/'^Pi^2|</.2p02 + 0(n2-5). 

2 

(64b) 

where 0{n^) is in the coarse-grained sense; it does not 
exclude the possibility of a correction term which may 
behave like [In (na^) + const.]. 

Like in the case of structureless bosons 2] , we express 
the particle number N in terms of the wave function, in 
order to solve it backwards for rj^- Up to the order n^-^, 
there are only two diagrams for TV = c\ci); they are 
the derivatives of Ti and Tia with respect to —fi. 



Solving the above equation for rja, substituting the re- 
sult back to the formulas we obtained so far, and solving 
dE/dN = fi, we get our results, summarized below. 



1 /2 

fi — E,n/2 + namn^ 1 + (32/3^7?) (n|af„) + h.o.c. 



1 + (128/15 V^)(7i|a^y^^ + h.o 



ai = r]a(f)i + ai 



■q^ = n]'^ 1 - (4/3^^F)(?l|af„)''' + h.o.c 



3 ^l/2 



(65a) 
(65b) 

(65c) 
(65d) 



a, 



TV, 



3/2 



l + (28/3V^)(nta^)'/' d'^/^ 



~ (128/15V^)(nTaL) ^ EAil<^2|'</'2 +0(r^^-'), 

2 

(65e) 

/3i234 = TVt [l + (16/3V^)(nTa^)'/'] </)'i234 + 0{n^n-^), 

(65f) 

/??2 = x^M2 + Pt^, E 'l'*i<l'*2Pt2 = 0, (65g) 

12 



x^ = -[l + _ ^q^2 + eq^) + h.o.c, (65h) 

= [(1 + 2a;q)(d'i02 + 0i4) + 9'i2] + h.o.c, (65i) 
7i'234;5 = 2v/]vfa;q0'i234(/'5 + h.o.c, (65j) 

where in Eqs. (j65f|) and ( |65j| ), the first four spin-momenta 
1234 must be an 1-cluster, thus excluding the cases in 



FIG. 7: Diagrams for the superfluid pairing function Fk 
(c-kCk) up to the order . 



which the sum of any two of them is ^ ^nam or 
zero. In Eqs. I|65h|) . (|65i() . and | |65j| ), the momentum 
q < min(l//, l/r„). 

As a consistency check, we calculate the sum of all 
the diagrams in Fig. |S1 The calculation is tedious but 
straightforward. The result is consistent with the value 
oiE- nN from Eqs. and (j65H|l . 



V. PHYSICAL OBSERVABLES OTHER THAN 
THE EQUATION OF STATE 

After a multistep process of logical deductions and an- 
alytical calculations, we have now determined the param- 
eters of the wave function [Eq. (QJ] to the leading order 
beyond mean-field, in the low-density regime. We are 
now in a position to determine various physical observ- 
ables, by evaluating their expectation values under the 
wave function. 

For every observable computed below, we first use 
Eq. (|29|l and the other information in Sec. lIV Bl to identify 
all the diagrams up to a certain order in the low-density 
expansion; we also note that there are no dead ends. 



A. Superfluid pairing function 

In this subsection we determine the superfluid pairing 
function, and discuss the range of densities in which the 
result is valid. 

We decompose it into a component parallel to the in- 
ternal wave function of an isolated molecule, and one 
orthogonal to it: 

Fk = (c_KCK> = fI + F^= r;j.0K + ^k, 

where Y.\^ '/'k-^K = 0- 

All the diagrams up to the order are shown in Fig.Q 
Retaining terms to the order n^, we get 

Fi =ai(l-|ai|2)- ^ (ajI?i+2«iCi^i)+0(^'-'). 

0<9<fcc 

(66) 

Using the established results of the parameters of the 
many-body wave function, we get 



1 - 



/ Q xl/2 , 

(n|a„j + h.o.c. 



30F 



n]'^^^, (67a) 
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28 



("T«m) 



3/2 



(d[/n-J2Pii\h\^h 



and the latter equation is simplified by Eq. (|18afl : 



28 



1/2 



OK'^). (67b) 



The accuracy in the parallel component of Fn is limited 
by that of the parallel component of aK ■ To increase the 
accuracy, one must determine E — fiN to higher orders. 

The amplitude of the component in parallel to 
(/)K is mapped directly to the well-known condensate- 
depletion formula 



{bo 



{ubal 



,1/2 



h.o.c. 



in the limit of point-like bosons. Here Ni, is the total 
number of bosons, and nt and at are the boson density 
and scattering length, respectively. 

The orthogonal deviation of -Fk from riF4>'K. is propor- 
tional to c^K - a function derived from the two-molecule 
zero-speed scattering wave function - at the leading and 
next-to-leading orders. The coefficients are also predicted 
exactly, in Eq. (|67b|l . The next step is to determine c?k 
from Eq. H18a|) . If the interaction between fermions is 
zero-range with a large and positive scattering length, 
we should substitute the four-fermion wave function cal- 
culated by Petrov et al into Eq. (|18a|l . to determine 
die. This step is elementary (involving no many-body 
physics at all) but tedious; it will be carried out later, by 
the author himself or by some other people. 

The significance of this deviation from the two-body 
function (/)k has been pointed out in the closing para- 
graphs of Sec. IIV Gl Quantitative comparison between 
our prediction and experiment is possible in the future. 

For ultracold fermionic atoms in two internal states 
near a wide Feshbach resonance (on the a > side), 
since our prediction is based upon the two-molecule scat- 
tering physics as studied by Petrov et al our pre- 
dictions are relevant to experiment in so far as Petrov 
et aPs results are relevant to experiment. There is al- 
ready much experimental evidence su ppo rting the value 
of am as predicted by Petrov et al most notably 

the cloud size of a trapped molecular Bose-Einstein con- 
densate. On the other hand, the two-molecule physics 
of Petrov et aVs is never exact except when the molec- 
ular density approaches zero; our low-density theory is 
in a similar situation. We are thus led to believe that 
the low-density results in this paper are applicable to ex- 
perimentally accessible densities, just like Petrov et aVs 
[ll| . So the low-density results in this paper are applica- 
ble when < kpa < 1. The more detailed upper bound 
depends on the particular experiment one wants to ad- 
dress, and on one's accuracy goal. Here kp = {Sir^ny/^, 
and a is the scattering length between fermions. 



I I 
I,-'' I 



I,'-' I 



FIG. 8; Diagrams for the superfluid four-fermion function up 
to the low-density order n^, when the sum of two external 
spin-momenta is ~ -y/n. The external points are distinguish- 
able. In the second diagram, the shaded circle stands for 
various internal structures (either N2 ~ 1, 2, or there is one 
bubble island that is not simplest). 



Having stated these, we must stress that our results are 
also valid beyond the range of validity of Petrov et al's 
\n\ in the sense that ours are valid for other interactions 
between fermions as well. For other interactions, the two- 
molecule scatte ring wave function will be different and 
^ 0.6a [H Il2l |. But the results in this paper still 
apply, provided that the conditions specified in Sec. Ill Al 
are satisfied. 



B. Superfluid four-fermion function 



Let Ki + K2 + K3 



K4 



(4) 



K1K2K3K4 



0, and 



(CK4CK3CK2CK1; 



Case 1: the sum of two of the four spin- momenta is 
zero, and the thermodynamically significant contribu- 



tion to F*^'') equals a product of two F*^ 



2)' 



Eq. ©. Here F^'^ = Fj 



K- 



s, according to 

(4) 

K,-K,K' 



-K' — 

FkFk' if K ^ ±K'. (If K = ±K', 4'Ik,k'.-k' = 0-) 

Case 2: for all i = 2,3,4, Ki and either have 
parallel spin, or have a nonzero total momenta ~ 0(n'^). 
There is only one diagram for F'^' up to the order n^'^, 
namely /3kiK2K3K4- 



5.(4) 

K1K2K3K. 



1 
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( 3 \l/2 



^KiK2K3K4 



(68) 

Case 3: the sum of two of the four spin- momenta is a 
small momenta ~ n^/'^. Consider for instance 

p(4)q ^ p(4) 

-^KiK2 — -^Ki+q/2,-Ki+q/2,K2-q/2,-K2-q/2' 

where qj ^Jn^a^ is a nonzero constant of order unity (in 
the low-density limit). If Ki 7^ ±K2, all the diagrams 
up to the low-density order are shown in Fig.|S| Using 
the part of the cluster-separation theorem that is already 
proved, together with the established result concerning 
7i234-5j easily show that the second diagram (con- 

taining the shaded circle), whose low-density order is 0.5, 
is proportional to (/)i02 up to the order n^. So 



^(,^)'»=(l-2|aip-2KP 

+ ??''</'i02 + O(ni-5f7-i) 



^12 
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where rf^ ~ -n}/'^. We first substitute the result of (3^2 
[Eq. (innj] to Eq. (|53b|l to determine A^^ accurately, and 
then evaluate the other terms to the leading order. At 
Ki = ±K2, there is appropriate delta- function singular- 
ity in ^"^1X2 ' accordance with case 1. Combining these 
subregions of the spin-momentum configuration space, 
we derive the final result 



(4)q 



12 



(4)q|| 



12 



F, 



(4)q_L 
12 ' 



p(4)q|| 
■^12 



h.o.c. 



,*^Wq-L 



9l(?>2, 



-^12 — "T 



1 -|- 2x 

^{di(j>2 + <l>id2) + gu 

l-a;2 



0, 



(69a) 



(69b) 



The physical meaning of these results will be clearer in 
coordinate space. 



C. Superfluid six-fermion function, momentum 
distribution, two-body reduced density matrix, etc 

The six-fermion correlation parameters 71*234.5 are an 
indispensible part of the consistent beyond-mean-field 
theory. This is not surprising, given the fact that 
three-boson correlation is needed to formulate the cor- 
rect beyond-mean-field theory of structureless bosons 
0- Leggett (and many others) notices that Bogoliubov 
wave function does not yield the correct beyond-mean- 
field equation of state of structureless bosons, since some 
"bare interactions" are not renormalized in favor of the 
bosonic scattering length, a problem he considers "a lit- 
tle disturbing" , which he then remedies by modifying 
the Bogoliubov wave function. 

Using the pseudopotential as the effective bosonic in- 
teraction, one can escape this problem at a price: the 
correct short- distance behavior of the bosons upon colli- 
sions is omitted. For the Fermi superfluid in the molec- 
ular Bose-Einstein condensated state, the price of this 
replacement is even higher: since all the four-fermion 
scattering physics is replaced by a molecular scattering 
length, one can no longer describe the many physical ob- 
servables that involve momentum scales much larger than 
-^nUrm such as the many-body effects on the fermionic 
momentum distribution and the superfluid pairing func- 
tion. 

In the approach of and this paper, such prob- 
lem is eliminated by strictly applying the power- 
counting formula to identify all the significant dia- 
grams. When this is done, one finds that six-fermion 
(or three-boson) correlation terms can not be altogether 
thrown away as people wish to do. In the same way 
as the last two subsections, one can derive the as- 
sociated nonvanishing superfluid six-fermion function: 

(cK4+q/4 • • • CKi+q/4C-K5-q/2CK5-q/2)- 

The momentum distribution of fermions, (c^Cp--), is 
another interesting quantity, for two reasons. Firstly, 



it is directly measurable. Secondly, it is predicted that 
there are two exact relations between the energy and the 
momentum distribution of the two-component Fermi gas 
with large scattering length '3, '3| , valid for all nonzero 
values of kpa (including the unitary regime), and also 
valid for finite temperatures, few-body systems, etc. Per- 
forming the analytic calculation of the momentum dis- 
tribution in the low-density regime, we can directly test 
these earlier predictions [3, H ■ The results obtained in 
this paper already enables an unambiguous determina- 
tion of the fermionic momentum distribution up to the 



rder 



,2.5 



It will be presented in a subsequent work. 



The two-body reduced density matrix is also observ- 
able. It remains to see if there is interesting physics in 
this quantity. It will be determined analytically in the 
low-density regime in a subsequent work. 

Another quantity to study is the fermionic Green func- 
tions ( "normal" and "anomalous" ) . They help us to un- 
derstand the dynamic behavior of the Fermi superfluid. 

The above list is still not exhaustive. 



VI. CONCLUDING REMARKS 

At low densities, there is a natural coexistence of the 
characteristics of a conventional Bose gas and that of a 
Fermi gas in the present theory. If the system is probed 
at length scales of the order ^ ~ l/^n^Om or longer 
length scales, one can hardly distinguish it from a usual 
Bose gas. The equation of state, the condensate deple- 
tion, the low-momentum part of the momentum distribu- 
tion of molecules {q ~ etc, all follow conventional 
Bogoliubov-Lee-Yang theory. If the observables that in- 
volve large momentum scale q ~ 1/rm are measured, 
however, one sees many features that are absent in a 
conventional Bose gas. The momentum distribution of 
fermions, the deviation of the superfluid pairing function 
from being proportional to the internal wave function of 
an isolated molecule, etc, are examples of such observ- 
ables. 

The Bose characteristics is not put in by hand, but 
forced on us by the fermionic wave function Eq. The 
ground state energy as derived in this paper is the abso- 
lute lower bound, as we have gone through great lengths, 
making all the possible adjustments of the wave function 
to minimize it. The wave function itself does not seem 
to be improvable any more, since all the possible correla- 
tions are included. It is the low density that makes higher 
correlations less important, but their effects are not auto- 
matically ignored - instead, using the power-counting for- 
mula Eq. H29|) . we have shown that most of them, except 
a few, are insignificant in the calculation of the leading 
order correction to the mean-field equation of state. 

Will such Bose characteristics be carried over to the 
high density regime? It seems that various neutral 
isotropic superfiuids (no matter BEC of structureless 
bosons, BEC of composite bosons, unitary Fermi su- 
perfluid, or BCS superfiuids) have very similar struc- 
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tures at length scales much longer than ^ (not neces- 
sarily cx rt^^/^). At shorter length scales, their differ- 
ences show up. However the similarity between the dilute 
BEC of structureless bosons and dilute BEC of composite 
bosons is kept at shorter length scales, until max(r,„,^), 
which is much shorter than ^. For a unitary Fermi super- 
fluid at zero temperature, there is only one length scale 
^ ~ l/kp ^ n^^/'^, so its similarity to other superfluids 
must be broken at a length scale somewhat longer than 
t 

To formulate the quantitatively well-controlled many- 
body theory for the Fermi superfluid in the unitary 
regime, we should probably take advantage of the above 
partial similarity, and do a momentum scale division at 
fee, where fee is roughly the highest momentum scale helow 
which the structure of the system is nearly indistinguish- 
able from other superfluids. The structure of the system 
at k > kc needs to be determined nonperturbatively, and 
smoothly matched to the low-momentum structure. 

Besides this project, we have of course the shorter run 
subject of extending the current low-density theory to 
higher orders beyond Lee- Yang. We need to first show 
that Wu's logarithmic term ^ is present in the equa- 
tion of state, and then determine the leading differences 
of the equations of state of molecular condensates and 
those of other BECs. More concretely, we need to deter- 
mine the value of the constant C as defined in Ref. 0. 
We derive confidence in these plans from two sources: 
the diagrammatic theory in this paper (the author has 
found diagrams in the expansion of E — which scale 
like n^^lnn), and Braaten et al's Effective Field Theory 
calculations 0|, which support the universality of the 
Wu term, at least for structureless bosons. 

In the same time, the fermionic momentum distribu- 
tion is of high interest. Here we have some exact theo- 
rems 0, 01 for the two-component Fermi gas with large 
scattering length, valid in the low-density as well as uni- 
tary regimes. They relate the momentum distribution to 
the energy, the pressure, and the rate of change of energy 
during a ramp of the fermionic scattering length [3, 
The confirmation of these theorems in a concrete low- 
density analytic calculation will push our understanding 
of such a novel Fermi system to a greater depth. This 
has been done in Ref. Q], for the dilute Base gas with 
large scattering length, although the theorems are only 
approximate in that Bose system 
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APPENDIX A: AN IDENTITY SATISFIED BY 



A(4) 



Proof of Eq. (gnj). Replacing Ki • • • K4 in Eq. (Pc|) 

by Ki + q/2, -Ki + q/2, K2 - q/2, and -K2 - q/2, 
respectively, using the Galilean symmetry of U, multi- 
plying the equation by '/>k2i summing over K2, noting 
the antisymmetry of (f> and t/)'-^-' , and noting that dumb 
momenta can be freely shifted, we get 



«2 

2 



T 



— 



^Ki+q/2 ^Ki+q/2 ^ '^^ 



,» ,(4) 
K+q.-Ka-q^KaKifKa+q/a'^^KKaKsKi 



and the infinitesimal e*^*^^ is omitted in comparison with 
(?V2 > 0. 

Using Eq. I)12d|l . we see that the first term of Eq. (|Aip 
approaches -(27ram/n) X]2l'^2p</>1 = -(47ra„i/r2)0Ki , 
when q ^ but q ^ 0. 

Obviously 1^4 -O14-P43 = ^i3- the second term of 
Eq. I|A1|I is rewritten as J^a ^^23 ^43'/'32^'''/'2- 

^2 -^2302 = Oj the third term vanishes. 

Since U and (j)* are smooth functions of momenta 
(when momentum is conserved by the subscripts), Tk = 



limn 



>0 



K3K, 



K -Ka.KjKi'pKa 



''KKoKsKi 



Obviously Tko- depends continuously on the momentum 
k, for a given spin a. So in the limit q — > 0, the last two 
terms of Eq. IjAip approaches T_Ki — Tki • 

Combining the above results, we see that when q — > 0, 
Eq. is obtained. 



APPENDIX B: DEFINITION OF SOME 
TECHNICAL TERMS 

The terms used in this paper follow convention as much 
as possible. However some terms in this paper have 
narrower meanings than convention, some are used in 
a broader sense, and some are new - due to the needs of 
the theory. Here we give the precise meanings of some of 
them (for use in this paper), in alphabetical order. 

Dead end: an internal part of a diagram that is con- 
nected to the remaining part via a single (de)generator. 
In the present formalism, dead ends are absent because 

it is zero 



ip) 

Ki-Kp- 



of a basic property of the coefficient a 
whenever the sum of a nontrivial subset of the p spin- 
momenta is zero. 

Diagram: an arbitrary collection of vertices and lines; 
each line has exactly two ends, and they must be attached 
to two vertices or the same vertex. The number of ver- 
tices is denoted by Nq, the number of lines by A^i, the 
number of independent loops by 7V2, and the number of 
disconnected parts - each part is connected, but any two 
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different parts are disconnected - hy N3. There is a sim- 
ple universal identity 

No- N1+N2- N3 = 0, for any diagram. (Bl) 

This identity is used in the derivation of two basic power- 
counting formulas, Eqs. (O and (|29|) . 

External point: a special vertex representing a 
fermionic operator (annihilation or creation) in a product 
whose expectation value we are computing. Each exter- 
nal point is attached by a single solid line associated with 
a certain spin-momentum. It is a small solid circle if it 
stands for a creation operator, or a hollow one for an 
annihilation operator. 

h.o.c: higher order correction (in the low-density ex- 
pansion) , with strictly greater value of low-density order 
R, in comparision with the term immediately preceding 
it. 

Low-density order i? of a quantity x: if the leading 
order approximation of x scales like n^OP in the low- 
density regime, the low-density order of x is said to be 
R. Symbolically we may write x ~ . If actually a; ~ 
n''^ln(na^) or things like that, we may still write x ~ 
in the slightly broader sense. Here n and r2 are the 
particle density and the system's volume, respectively, 
and O is assumed to be large enough to avoid finite-size 
effects. A higher low-density order means a greater value 
of R, and usually indicates a smaller magnitude in the 
low-density regime. 

Nontrivial subset of M items: a nonempty proper sub- 
set, namely with Mi elements such that 1 < Mi < M—1. 

O(n^): abbreviation of Oin^n°). 

0{n^n''^): a symbol describing how a quantity (in the 
low-density regime) scales with n and fl. The (often di- 
mensionful) coefficient roughly independent of n and Q is 
not shown in this symbol. Things like ln{na'^)Q'^ are 
sometimes possible, and not excluded by this symbol. 

Thermodynamic order Q of a quantity x: if x scales 
like in the thermodynamic limit (in which n is held 
constant as usual) , the thermodynamic order of x is said 
to be Q. Symbolically we may write x ~ fi*^. This notion 
is obviously orthogonal to the low-density order; it is also 
valid beyond the low-density regime (in contrast with the 
low-density order). 

No contradiction between x ^ il'^ and x ~ exists 
- there is no contradition between "the height of this 
rectangle is Q" and "the width of the rectangle is R" . 
In fact we may use a plane, and map a; to a point with 
coordinates (i?, Q) on the plane. 



APPENDIX C: PARTIAL PROOF OF THE FIRST 
CLUSTER-SEPARATION THEOREM 

First consider _Fk = (c_kCk)- The theorem states 
that there exists a coefficient rj = 77g gn°-^ -I- rj^ qTi^, such 
that i^K = v4'K.i with an absolute error of the low-density 
order n^'^. The proof follows. 



K2 

(a) o^^^^^^o (b) o-C^.-o (c) o-^/\ 
K4 

FIG. 9: Diagrams for i ^ t^KiK2K3K4(cK,CK4CK3C„Ki>- 
K2K3K4 

(a): the total diagram. The shaded circle stands for the 
ground state expectation value {cl^^CiQ^Cj^^c_ji^) . (b) and 
(c): two individual diagrams, whose low-density orders are 
1.5 and 2, respectively. In (c), if the whole island is reduced 
to a point, the dotted line connects this same point and forms 
a loop. 

Since the ground state expectation value of 

H = H - fiN (CI) 

is stationary under any infinitesimal change of the many- 
body state (the number of particles may also be changed 
slightly), we can easily show an equation of motion 

{OH - HO) = (C2) 

for any product of fermionic operators, O. Substituting 
Eq. (121 and O — c_KiCki into this formula, we get 

5x2^2 = i ^ C^1234(4c4C3Ci) + i ^ C/2134 (4C4C3C1) , 
2 234 234 

(C3) 

-DK1K2 = {kf — Em — 2^)(5kiK2 + 2^Ki,-Ki,K2,-K2- 

(C4) 

where fi = fj, — Ejn/2 is at least of the order nP-^ 

The first term on the right side of Eq. (|C3|I has a di- 
agrammatic representation in Fig. ^ In Fig. EJa) the 
total diagram is shown; the big shaded circle stands for 
the expectation value of a product of fermionic opera- 
tors, which in the present case is (cKjCj^^Cj^gC.K:). In 
Fig-IHb)(c) we show, for illustration, two associated in- 
dividual diagrams, whose low-density orders are 1.5 and 
2, respectively [Eq. note that iV;^ = 1 for Fig.EIc), 

since its skeleton diagram is a point and a dotted line 
whose both ends are attached to such a point. 

We now show that the low-density orders of the 
infinitely many individual diagrams associated with 
Fig. Ela) have a lower bound, 1.5. In Fig. Ufa), the 
piece on the left of the shaded circle must be in the 
same island as the solid line on the right, because the 
dotted lines (if there are any) can not carry spin or any 
large momentum, and the separation of the two pieces in 
two different islands would contradict spin-momentum 
conservation. Thus we have a single external island e. 
One spin-momentum (ie K2) enters the shaded circle, 
three leave it, and the associated four fermion opera- 
tors are normally ordered, so K2 must be attached to 
a degenerator in the big shaded circle; so Pi ■* > 2 and 
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Ki+q/2 



■K2-q/2 Ki+q/2=-<^ K2-q/2 



-Ki +q/2 o^W^o -K2-q/2 -Ki +q/2 \ -K2-q/2 

FIG. 10: Some residual terms in the equation of motion of 

(c-Ki+q/2CKi+q/2C-K2-q/2CK2-q/2) [see Eq. lEU]. In the 
second diagram, "f marks a later "time" than the other legs 
of the shaded circle, as the six fermion operators are normally 
ordered. The low-density orders of these terms are at least 1. 



ii"' - Ll'^' = Pi ' + 2>i. So Pe> 6. 
Note also that — 1. So we deduce from Eq. (|31|l that 
R > 1.5. 

The last term of Eq. (jC3p is similar. So the func- 
tion {DF)k = X^K' -Dkk'Pk' is at least of the order 
n^'^. The spectrum of D has a gap equal to \Em\', the 
eigenvalue associated with 0k is — 2pt, and all the other 
eigenvalues are at least \Ejn\ — 2pt > \Em\/2 ~ n° (see 
Sec, nm for our basic assumptions about U) . Now make 
an orthogonal decomposition: Fk — ?70k + (where 
Sk'^K'^k = 0). Since D is hermitian and (^k is an 
eigenvector, {DF)k = —2'jlri(j)-i^ + {DF')k is also a sum 
of two orthogonal components, and the low-density order 
of each of them is at least 1.5. Since the leading order 
contribution to Fk is just an, the low-density order of 
r) is 0.5 (Sec. IIVC(I . So F^ is at least of the order 
and (as a byproduct) 



,1.5 



^ = - Em/2 

which is needed in Sec. IIV Fl 



(C5) 



Now turn to PkiKs 



*'^-Ki+q/2'^Ki+q/2'^-K2-q/2 



'^K2-q/2)' ^ function of K1K2 (where neither Ki + 
K2 nor Ki — K2 is zero or ~ \/n), for a given nonzero 
momentum q (9 ~ V^)- The low-density order of this F 
is 0. 

The cluster-separation theorem states that there exists 
a coefficient r]q such that PkiK2 = '?q'/'Ki0K2 +-P"kiK2' 
where the low-density order of F' is 1 or higher. To prove 

it, we substitute O = C_Ki+q/2CKi+q/2C-K2-q/2CK2-q/2 



into Eq. I|C2|1 and get 



(22)q 



■ F, 



q 



4 

1=1 



4 

i=l 



where -ffS^^.-KiK^ 



(2)q 
KK' 



En 



2p^-f gV4)J] 



^ A ff(2)q 
, + '^KiKi^K2Ki' 

KK' + 5C/k,-K,K',-K', 



and Ci and are eight residual terms whose low-density 
orders are at least 1 (see below). In Fig. EH are shown 
the diagrams for ei and /i. The other e^'s differ from 
ei in the two external points that are linked to the solid 
side of the diamond; but when Ki + q/2 and — Ki -I- q/2 
(or K2 — q/2 and — K2 — q/2) are linked to the solid side 
of the diamond, we get some terms on the left side of 



Eq. ljC6|l . The other /^'s differ from /i in the external 
point that is linked to the solid side of the diamond. 

Now apply Eq. H31|l . For ei, all the external points 
must be in the same island, so /cxt = 1; also Pcxt > 4; 



so i? > 1. For /i, /ext < 2, Pit' > 2, P^.V > 6, and 



,(+) 



(+) 



Pcxt — Pcxt 



p'^^ > 



soi?>l. 

The spectrum of iJ^^^i has a gap: the lowest eigenvalue 
is of the order and is associated with the eigenvector 
0K, and all the other eigenvalues are at least \Em \ — 2'jl + 
> \Em\l2. The two terms of if(22)q (^^^ ^(2)q) 
independently on Ki and K2, respectively. So i/(22)q 
has a single eigenvalue of the order v}, associated with 
the eigenvector 0Ki 0K2 ! all the other eigenvalues are at 
least \Ejn\ — 4pt-|- /2 > \Em\/2. So the components of 
^KiK2 '^hich are orthogonal to (/'Ki'^K2 must be of the 
order , in order for i/(22)q^ j^g q£ ^^le order . 



On the left side of Eq. lfC6|) . K'^ (or K2) may take val- 
ues close or equal to ±K2 (or ±Ki). But the associated 
contributions are not important. When K.'i =F K2 ~ y/n, 
the associated volume in the K'j^-space is reduced by a 
factor of the order n^-^. When K'j^ =F K2 is exactly zero, 
the thermodynamic order of Pk'jK2 increases by 1, but 
this is accompanied by the loss of an independent dumb 
momentum; meanwhile the low-density order of Pk'j^K2 
becomes 1, so the net effect on Eq. (IC6|) is of the same 
order as those of and fi. 

For more 1-clusters {M > 3), each of which contains 
two spin-momenta, we still find that the equation of mo- 
tion is dominated by the quasi-independent dynamics of 
individual clusters, with small "driving terms" [like the 
ones on the right side of Eq. (|C3|) or Eq. (|C6I) ]. so that F 
(at given small momenta q; ~ ^/n) is still proportional 
to a product of single-molecule internal wave functions, 
with relative error of the order n^. 



APPENDIX D: PARTIAL PROOF OF THE 
SECOND CLUSTER-SEPARATION THEOREM 

Consider any diagram V in the low-density expansion 
of Pk = (c-kCk)- Since there are only two external 
points [L'^^^ — 2, = 0], there is only one external 



island, and /cxt 
deduce that R = 



P 



(-) 



= 1, P^Vj^ = P^J + 2. From Eq. (EU we 

1/2 + Pi-V2 + (Pint/4 - /int) + A^^/2. If 

0, the external island is of a unique structure - a 
single small solid circle linked to the two external points 
via two solid lines; moreover, this island must be the 
whole diagram, since dead ends are absent. If P^~^ > 2, 
P > 1.5. So Pk = «k plus correction terms of the order 



n 



1.5 



But we proved in Sec. IIVDI that there exists an 77 
such that Pk = ?y0K plus correction terms of the order 
n^-^. So aK = '70k plus correction terms of the order 
The simplest case of the second cluster-separation 



n 



1.5 



theorem is thus proved. 
Now turn to 

/5kiK2 = /^Ki+q/2,-Ki+q/2,K2-q/2,-K2-q/2 (Dl) 
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FIG. 11: First row: diagrams for F^^j^^ up to the order 
nP'^. Second row: diagrams for J'kiK2K3 order 
n*^. Each dashed line represents a geometric series associated 
with a variable number of simplest bubbles. The first diagram 
corresponds to the first row of Fig. El 



as a function of K1K2 (where neither Ki + K2 nor 
Ki — K2 is zero or ~ \Ai)j fo'" fixed small q{q ^ V^)- We 

First 



extract its property from -FkiK, i^^*^ IIVD|I 
identify all the diagrams in the low-density expansion of 
F satisfying R = or R — 0.5; using Eq. we find 

that in such a diagram, there are two external islands 
(each of which has Pi — 2), the internal islands are sim- 
plest bubbles, and TV^^ = for i? = 0, iV^ = 1 for R = 0.5. 
These two sets of diagrams are shown in Fig. ^2 (first 
row). Some diagrams have the same structure after the 
simplest bubbles are suppressed, and their sum can be 
represented by a single diagram; each diagram in Fig. 1111 
is such a sum; a dashed line represents a geometric series 
associated with a variable number of simplest bubbles. 
See the first row (^kk') '^^ Fig. for the content of the 
first diagram in Fig. 1111 



^KiK2 ^ "^KiKa 



(D2) 



^KiK2 



/^KiKa + X] '^KiKa q^K^Ks o'^KiKa 
K3K4 



(D3) 

where 1/2 is the symmetry factor of a simplest bubble, 
and the small momenta flowing through two adjacent 
bubbles have opposite signs (and directions). In the sum- 
mation over an internal spin- momentum, eg K3, K3 may 
be close or equal to ±Ki or ±K4, but the corresponding 
contributions are negligible for our purpose, since their 
associated volume in the dumb-momentum space is re- 
duced by a factor of the order n^'^. Now introduce ma- 
trices /3q and ^q, whose matrix elements are defined by 
Eqs. Hm|l and HU3() . respectively: (/3q)KiK2 = /^KiKs 
and similarly for Aq. So 



An 



where 



+ /3q/3^/3q/4 + /3q/3^/3q/3^/3q/16 + • • • 



Pl(3JA (D5) 



is a hermitian matrix, and all its eigenvalues are nonneg- 
ative. Also the geometric series in Eq. (ID4II should be 



convergent, so the eigenvalues of should be less than 1 



Since A^Aq = (l-rq)-i/?t/J^(l-rq)-i =4rq(l-rq)-^ 
1 + AjjAq = [(1 + rq)(l - rq)"i]2. But all the eigen- 
values of the hermitian matrix (1 -I- rq)(l — fq)^^ are 
greater than or equal to 1, according to the above prop- 

-1 



erties of rq. So (1 + 7-q)(l — r^Y 



1 



AqAq. 



[The 



square root of any positive definite hermitian matrix h, 
VTi, is defined as a hermitian matrix with the same set 
of eigenvectors as h, but its eigenvalues are the (posi- 
tive) square roots of the corresponding eigenvalues of h; 
this definition is unambiguous even when the spectrum 
of h has degeneracy; y/h depends smoothly on h.] So 

1 



1 



1 



2(1 -rq)- 



2A„ 1 



and 



,(> 



AqAq 



(D6) 



So /?q is a smooth matrix function of the matrix A^ 
K1K2 



is the i? = 0.5 correction to F (the second 
diagram in the first row of Fig. 111(1 . 
??q0Ki0K2 +0{n/n). So 



and i^^^K2 



A"^ 

^KiK2 



77q0Ki(/'K2 +0(nO■•V^^). 



(D7) 



Substituting this into Eq. (|D6|I . we get [after using 
Eq. pUct ] 



^KiK2 ^ 



2?7q 



l + v/l + 4|r7q|^ 



»Ki<PK2 



+ 0{n"-^n-'). (D8) 



This is still one step from the desired result, namely 
/^KiKa ~ '7q0Ki0K2 +0(n^f2~^) for some r^q. So we need 
to consider the "1-loop" term -BkiK2' and show that it is 
factorizable at the leading order (~ n°-^). But this term 
contains the dispersed vertex vm (or fjij^). So we need 
to first show that wm, and by the way wim, wnm, etc, 
are factorizable at the leading order. 

(C-Ki+qi/2CKi-|-qi/2 



So we turn to F^^'^jl^^^ = 



-0.5 



C-K2-hq2/2CK2-Hq2/2C-K3-|-q3/2CK3+q3/2; ~ n (whcrC 

qi + q2 + qa = 0, and gi,2,3 ~ V^)- Using Eq. we 
see again that diagrams at the lowest two orders have the 
same external-islands structure. The leading expression 
is the first diagram in the second row of Fig. ^2 It has 
three dashed lines, linked to the external islands associ- 
ated with qiKi, q2K2, and qsKs, respectively; suppose 
that the numbers of simplest bubbles on them are si, S2, 
and S3, respectively. Because the small circles contained 
by a dispersed vertex must be of the same color, si_2,3 
must be all even (for vm) or all odd (for fjn). After 
straightforward derivation, we get 

pqiq2q3 \ ^ /^qi /^q2 /^qa ^qiq2q3 

^KiKaKa ~ ^KiKJ '-^K2K^ "-^KsK;, 'K'^K^K^ 

■^1,2,3 

I 1 /iqi ^*qiq2q3 

^ o ^KiK'j"^K2K;j^K3Kf, 'K'jK^K^ 



K', 



(D9) 
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-^^^^ /KiKaK;,— 'Ki+qi/2,-Ki+qi/2,...,K3+q3/2 -K3+q3/2' 

the matrix Gq = (1 — /3q/3^j/4)^^, G^^j^^ is its matrix 
element, = — q,, and B^^"^^^ ~ n° is the next-to- 
leading order contribution (the last three diagrams in the 
second row of Fig. [TT|) . 

The first two terms on the right side of Eq. (|D9|1 are 
the leading-order expression of F, but the second term 
is proportional to 0Ki0K2'/'K3 at the leading order, ac- 
cording to Eq. IID7I) : since F has this same property, we 
deduce that the first term on the right side of Eq. ljD9|l 
is also proportional to </'Ki0K2 0K3 at the leading order. 

The linear kernal in the first term is the direct product 
of three matrices: JC = Gq^ (g) Gq^ <^ Gq^ . Using Eq. (|D8p , 
we find that at the leading-order, {Gq4>)K — •^q'/'Ki where 
Aq = 2/(1 + ^1 + 4|ryqp) > 0, and aU the vectors or- 
thogonal to are unchanged under the action of Gq. 
So the hermitian matrix K. is positive definite, and has 
0Ki0K2'/'K3 as one of its eigenvectors at the leading or- 
der. 

Combining the two findings in the above two para- 
graphs, we deduce that TkiKs'Ics ™ust be proportional 
to </'Ki0K2 0K3 at the leading order. 

Now turn to ctj^^'ii^K^Ki ' namely the coefficient of the 
dispersed vertex wim. There are two diagrams, I?i and 
I?2, for J^KiKa^KsKi ^^1^ leading order. The skeleton of 
each of them is a tree (iVj — 0). One of them, Vi, con- 
tains two Will's (or two fill's, or vm and wjn) so it is 
factorizable at the leading-order according to the estab- 
lished properties of vu and fm. The other diagram, I?2, 
contains a single wim (or wJin) together with variable 
numbers of vu's and wJi's. Because both Vi and F are 
factorizable at the leading order, I?2 is also factorizable. 
Using the same logic as in the case of vm above, we 
can then deduce that fmi is factorizable at the leading 
order. 

Continuing this analysis, we see that for any M <C 
-^i Q^k^^'km (^^^ coefficient of the dispersed vertex 
''^l . . . l) is proportional to n^£]^(/)K„ at the leading order, 

M 

for any given set of small momenta q's. 

So we can now start the second round of logic, to 



strengthen the above result. 

First consider vu. In Eq. I|D2|1 . the term i?KiK, corre- 
sponds to the "1-loop" {N2 = 1) diagram in the first row 
of Fig. 111! and must be factorizable (ie, proportional to 
0Ki 4''K2 at any given small q ~ ^/rl) at its own leading 
order, according to the estabhshed properties of Vu and 
Will. But B is just the next-to-leading order contribution 
to F. We have proved in Sec. IIV 1)1 that F is factorizable 
at the lowest two orders. So the term ^KiKa ™ust also 
be factorizable at the lowest two orders, namely there 
exists an 77^, such that ^KiK2 ~ ^ /g^K i^Ka +0{n^il~^). 
Substituting this result to Eq. (|D6(I . we see that for 

< = 2</(i+yrT^^), 

Pk.k, = Vq4>K,cf>K, + 0{n'n-'). (DIO) 
We then turn to wm. In Eq. IID9|) . the term B cor- 
responds to the N2 — 1 diagrams in the second row of 
Fig. ^2 which consist of the dispersed vertices we have 
studied. B is thus factorizable at its own leading order. 
But F is factorizable at the lowest two orders, and B is 
just the next-to-leading order contribution to F. So the 
sum of the first two terms on the right side of Eq. ljD9p 
is factorizable at the lowest two orders. But the second 
term itself is factorizable at the lowest two orders, accord- 
ing to the above estabhshed property of ^kk'- the 
first term is factorizable at the lowest two orders. Note 
finally that Gq is determined by /3q, so (Gq0)K = •^q'/'K 
with a relative error of the order n^, for some A^. So we 
find that 

7SX^3 = <,q2q3'^K,0K2 0K3 (DU) 

with a relative error of the order n^, for some '/qjqjqg- 

Extension to more clusters (each of which contains two 
spin-momenta) is trivial. In general, we have thus proved 
that 

«l^f.'^^™^=<....q.>K,---</'K.. (D12) 

with a relative error of the order in the low-density 
limit, for any 1 < Af ^ N, provided that ± Kp are 
not zero or ~ y/n (for all 1 < < p < M). 
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